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ABSTRACT 


The Spitzer Survey of Stellar Structure in Galaxies (S^G, Sheth et al. 2010) is a 
deep 3.6 and 4.5 /im imaging survey of 2352 nearby (< 40 Mpc) galaxies. We describe 
the S^G data analysis pipeline 4, which is dedicated to 2-dimensional structural surface 
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brightness decompositions of 3.6 /im images, using GALFIT3.0 (Peng et al. 2010). Be¬ 
sides automatic 1-component Sersic fits, and 2-component Sersic bulge -h exponential 
disk fits, we present human supervised multi-component decompositions, which include, 
when judged appropriate, a central point source, bulge, disk, and bar components. Com¬ 
parison of the fitted parameters indicates that multi-component models are needed to 
obtain reliable estimates for the bulge Sersic index and bulge-to-total light ratio (S/T), 
confirming earlier results ( Lanrikainen et aL]|2007 ; Gadotti|[2008 ; Weinzirl et al. |[^009D . 
In this first paper, we describe the preparations of input data done for decompositions, 
give examples of our decomposition strategy, and describe the data products released via 
IRSA and via our web page (www.oulu.fi/astronomy/S4G_PIPELINE4/MAIN). These 
products include all the input data and decomposition files in electronic form, making 
it easy to extend the decompositions to suit specific science purposes. We also pro¬ 
vide our IDL-based visualization tools (GALFIDL) developed for displaying/running 
GALFIT-decompositions, as well as our mask editing procedure (MASK_EDIT) used in 
data preparation. In the second paper we will present a detailed analysis of the bulge, 
disk, and bar parameter derived from multi-component decompositions. 


Subject headings: galaxies: spiral —galaxies: kinematics and dynamics —galaxies: 
structure 


1. Introduction 


How and when did the baryonic mass assemble into galactic disks? How does the fraction 
of mass confined into bulges evolve over time? How common are galaxies that have no classical 
bulges, i.e. bulges that have their origin in the early mergers of dark matter halos and baryonic 
disk systems? These are difficult questions to answer because galaxy evolution involves secular 
processes such as gas accretion via filaments, where mass presumably ends up in bulges or disks, 
or internal dynamical evolution, such as the formation of bars which further re-distribute matter 
in galaxies. Galaxies in the local Universe are the present day manifestations of this evolution and 
hence provide important clues on the evolutionary processes which took place in the past. 


The Spitzer Survey of Stellar Structure in Galaxies (S^G, [Sheth et al. (2010)) provides an 
excellent data base with which to measure the stellar mass distribution of galaxies in the local 
Universe. It is a survey of 2352 galaxies observed in the mid-IR at 3.6 and 4.5 /im, wavelengths 
that are largely unaffected by internal extinction (Draine &; Lee 1984), and trace mainly the old 
stellar population (Pahre et al. ~| |2004[ however see also |Meidt et al.||2012] and [Driver et al. ||2013| ) , 
so that the mass-to-luminosity (M/L) ratio in these bands is nearly constant inside the galaxies 


(Peletier et al. 2012). This is particularly important for deriving the properties of bulges and 


disks, because dust and star formation are more pronounced in the disks than in the bulges, which 
in the optical region affect their relative M/L-ratio and thus the relative fraction of the bulge light 
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(Driver et al. 2013). Dust and star formation are significant also in the bulges of late-type galaxies 
dFisher |2006l |Gadotti|200Tt . The S^G images are deep, reaching azimuthally averaged stellar mass 
surface densities of ^ 1 Mq pc“^, where the baryonic mass budget at least in spiral and irregular 
galaxies is typically dominated by atomic gas. S^G covers a large range of galaxy magnitudes (over 
three decades in stellar mass), which makes possible to study both late-type dwarfs and bright 
galaxies in a uniform manner, and to study when the disk instabilities such as bar formation start 
to play an important role. Our sample extends to lower galaxy luminosities than most previous 


samples in which bars have been studied (Barazza et al. 2008; Sheth 2008; Nair et al. 2010 


Melvin et al.||2014 ). Besides galaxy mass, another central factor affecting its structural evolution is 
its environment ( van der Wel||200'8 ; Kormendy fc Bender||201^ ; Weinzirl et al.|[2014 ). S^G includes 
galaxies up to 40 Mpc and covers a wide range of different galaxy environments, including several 


galaxy groups and the Virgo and Fornax clusters (see Fig. 2 in Sheth et al. 2010). 

Plenty of information for the S^G sample is already publicly available via the IRSA archive. 
The data have been processed through Pipeline 1 ( Regan et al. |[20T4 ) (hereafter PI) which makes 
mosaics of the observed individual frames. Pipeline 2 (Muhoz-Mateos et al. P2) which makes 

masks of the foreground stars and image defects, and Pipeline 3 (Muhoz-Mateos et al. IIMII P3) 
which measures the basic photometric parameters like the galaxy magnitudes and concentration 
indices. In Pipeline 4 (P4), described in this study, we decompose the two-dimensional ffux distribu¬ 
tions of the images into several structural components using GALFIT (Peng et al. 2010). Because 
even the mid-IR wavelengths are not completely free of such contaminants as hot dust, mass maps 
are also created for the images in Pipeline 5 (P5, Querejeta et al. 2014). The galaxies in S^G have 
been visually classified at 3.6 /im by (Buta et al. [2010 , 2014), and we use these classifications in 
the present study. Optical images are also available for the majority of the S^G sample (Knapen 


et al. 2014). 


For all of the S^G galaxies for which the image quality is good enough (e.g. no superposed 
bright stars, or image defects), we provide 1-component single Sersic, 2-component bulge-disk 
(Sersic -h exponential), and multi-component decompositions, fitting up to four separate structure 
components. Our main goal is to estimate the parameters of the bulge and the disk in a robust 
manner, which is the motivation for our decomposition approach. In particular, it is important 
to include bar-components in the decompositions because the ffux of the bar is easily mixed with 
the ffux of the bulge (Laurikainen et al. 2006). Our bulge is defined as a ’photometric bulge’. 


including the ffux in excess of that in disk and bar components; the decompositions themselves 
do not make assumptions about the physical nature of the bulge, whether a rotation supported 


classical bulge or a disk star formation/bar vertical buckling related pseudo bulge (see Kormendy &; 
Kennicntt||2004 ; Athanassonla |[200^ . To measure the scale lengths and central surface brightness 
of the disks in a uniform fashion, an exponential function is used whenever possible, instead of a 
generalized Sersic function. It is well known that galactic disks can have more than one exponential 
sub-section (Freeman 1970; Erwin, Beckman & Pohlen |2005[ ). In this study we handle this in 
a fairly conservative manner: two separate functions (added together) are used to fit the disk in 









































































-4- 


galaxies where distinct inner and outer components of different surface brightness are present, but 
not in all cases in which a disk break (’truncation’ or ’anti-truncation’) of some degree has been 
reported in the literature. Our multi-component approach is similar to those used previously by 


Laurikainen et al. (2005, 2007, 2010), Gadotti (2009), and Weinzirl et al. (2009). Our motivation 


for offering also the single Sersic and bulge-disk decompositions is that they are routinely used in 


large galaxy surveys and high-redshift studies (HauBler ||2013; Lackner & Gunn ||2012; 

Cameron et 

al. 2009; 

Allen et al. 

2006; Driver et al. 

2006, 

2013). Although single Sersic fits are not good 


tracers of the properties of bulges, they are still useful in gross classification of galaxies. 

The decomposition results, released via IRSA and our web-page, are given in such a manner 
that they can be easily extended having different scientific goals in mind. The decompositions were 
done via GALFIDL, which consists of IDL-based tools for displaying and running GALFIT (see 


Sect 2.4). It is important to note that due to the large amount of work involved, P4 was started 


as soon as the first PI data was available. Because of this we did our own mask editing, and 
orientation and sky background estimatioi|^ These masks form part of the final P2 masks. Due 
to later changes in PI, part of the images used in P4 contain minor shifts (or differ in size by 1-2 
pixels) compared to the finalized PI images in IRSA. Rather than repeating the time consuming 
GALFIT decompositions with the updated images, we provide together with the decomposition 
output files the sky subtracted data and mask images we used. 


In this paper, we describe the decomposition method and model components, the preparation 
of the data for decompositions, and concentrate on illustrating our philosophy behind the construc¬ 
tion of the final multi-component decompositions. The results published in tabular form include 
the outer disk orientation estimates, Sersic parameters from the 1-component fits, and the final 
parameters from multi-component decompositions, together with a quality flag for each galaxy. 
The data products released via IRSA include the GALFIT output files, and all the input fits-files 
needed for repeating and refining the decompositions. The P4 web pages illustrate the same models 
in pictorial form, and also provide the GALFIDL code and documentation. (The IRSA products 
and the P4 web page are described in the two Appendixes). Analysis of the derived bulge, disk, 
and bar parameters will be presented in paper 2 (Salo et al, in prep.). 


2. Decomposition Pipeline 
2.1. Decomposition method and model functions 


Our decompositions use the GALFIT-software ( Pengetal. [2002 , 2010), which has become the 
de facto standard for detailed two-dimensional structural decompositions. It relies on parametric 
fitting, using the Levenberg-Marquadt algorithm to minimize the weighted residual between 


^The derived sky background values and orientation parameters turned out to be in very good agreement with 


P3, see Sect 2.2.3 
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observed (OBS) and model (MODEL) images, 


X 


1 [OBS{x,y) - MODEL{x,y)f 


( 1 ) 


The sum is taken over all the used (non-masked) image pixels, and cr(x, y) indicates the statistical 
uncertainty of each pixel (sigma-image). The model image consists of a sum of model components, 
i.e for bulge, disk, bar etc, convolved with the image Point-Spread Function (PSF-image). Note 
that the reduced is used, with N denoting the degree of freedom, equal to the number of fitted 
pixels minus the number of free parameters in the fit. 


GALFIT is extremely versatile in its selection of model components. Basically the user defines 
for each component its ’radial’ profile function, giving the surface brightness '^(r) at each isophotal 
radial coordinate r. The isophotal coordinates are most commonly defined in terms of generalized 
ellipses (Athanassoula et al. ][i99§, 


r{x\y') = 


\x' — Xq\^~^‘^ + 


y -yo 


y 


C+2^ 


C+2 


( 2 ) 


Here xo,2/o defines the center of the ellipse, q = h/a is the ratio between minor and major axis 
lengths. The x'^y' denote coordinates in a system aligned with the ellipse, with the major axis 
pointing at the position angle PA. For pure ellipses C = 0, while C > 0 indicates boxy and C < 0 
disky isophote^ For the pipeline decompositions, simple elliptical isophotes (7 = 0 are used for all 
components. Besides generalized ellipses, GALFIT provides several alternatives, such as definition 
of isophotal shape via azimuthal or bending modes, or via coordinate rotations, which would form 
a natural basis for detailed modeling of e.g. logarithmic spirals. To keep our models relatively 
simple (and uniform over the wide range of angular sizes and surface brightnesses spanned by the 
sample), we have not used these advanced GALFIT features. Keeping the models simple makes the 
interpretation of the observation minus model residuals more straightforward (see the NGC 1097 
examples in Sheth et al. ( 2010D ). 


The pipeline decompositions use five different choices for the model components/radial func¬ 
tions: 


1) The bulge component is described with a Sersic profile ( ^^sersic’ 


E(r) = Seexp (-k - l]) , 


( 3 ) 


^Note that in the original notation of 


Athanassoula et ah 


(1990) the exponent ’C+2’ was denoted with c, a pure 


ellipse thus corresponding to c = 2. Similar notation was used also in e.g. Gadotti (2011). However, we will here 
follow the notation of GALFIT ( Peng et al. |2002 2010) 
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where Eg is the surface brightness at the effective radius Re (isophotal radius encompassing half of 
the total ffux of the component). The Sersic-index n describes the shape of the radial profile, which 
becomes steeper with increasing n. In particular, n = 1 corresponds to an exponential profile and 
n = 4 to a de Vaucouleurs profile. The factor is a normalization constant determined by n. In 
GALFIT the corresponding -function is used, with the integrated magnitude mijuige a 

free parameter (instead of Eg). 

2) In the case of low or moderate inclination, the disk component is described with an infinites¬ 
imally thin exponential disk ( ‘^expdisk^^) , 

E(r) = T.oq~^ exp(-r//i,^), (4) 

where Eg is the central surface brightness of the disk observed from the perpendicular direction and 
hr denotes the exponential scale length. In this case the q — cosi, where i is the disk inclination. 
Assuming no extinction, is the projected surface brightness at the sky plane. The ‘^expdisk^’- 

function in GALFIT is used, with integrated rridisk = —2.5 log]^Q(27ri;o hr‘^) as a free parameter 
(instead of Sq). Note that in cases that had more than one disk component, the inner disk was 
sometimes fit with a sersic or ferrer2 function, to allow the profile to drop faster than with expdisk. 


3) For a nearly edge-on disk (apparent axial ratio q < 0.2), the function ( ‘‘edgedis¥^) 


^{rx,rz) = [j^ ) sech^{rz/hz), 


( 5 ) 


is adopted, where and are the (positive) distances along and perpendicular to the apparent 
major axis of the disk, and Ki stands for a modified Bessel function. This function corresponds 
to the line-of-sight (viewing along the disk plane) integrated surface brightness of a 3D luminosity 
density distribution ( [van der Kruit fc Searle |[l981 ) 


A(r^,r^) 


2hr 


exp{—rx/hr) sech^{rz/hz). 


( 6 ) 


4) For a bar component a modified Ferrers profile ( ‘Jerrer2^^) is assumed. 


S(r) 


So [l - [rlvoutf r < Tout 

0 r ^ T"out 


( 7 ) 


Here Tout defines the outer cut of the profile, while a defines the sharpness of this cut. The parameter 
/3 defines the central slope of the profile, and Sq is the central surface brightness (in the plane of 
the sky). 


5) When the galaxy contains an unresolved central component it is fit with a PSF-convolved 
point source ( ^^psf’^). In this case the free parameter is the total magnitude rupsf. Typically this 
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component, if present, is not an active or starburst nucleus, but rather a small bulge with angular 
size so small that it cannot be resolved in the S^G images {Re < FWHM = 2.1" of S^G images). 


For the decomposition pipeline we chose to do three types of decompositions: 1) one-component 
Sersic-fits, 2) two-component bulge-disk decompositions using Sersic-bulges and exponential disks 
(or edge-on disk if appropriate), and 3) multi-component ’final’ decompositions, optionally with 
additional bar, disk and central components (the level of complexity of the models is discussed in 
more detail in Section |^. The first two types of models are made in an automatic manner, while 
the final models always include human judgment about what components should be included. 


2.2. Preparation of data for decompositions 


2.2.1. What is needed? 


The S^G data analysis Pipeline PI (Regan et al. 2014) provides image mosaics in both 3.6 


and 4.5 /im, accompanied with weight-images^ which indicate for each pixel location the number of 
original frames covering it. Together with the header information, these weight images provide the 
means for producing the sigma-images used in GALFIT. 

Before decompositions can be started, frames masking the foreground/background objects and 
various image defects are needed. Additionally we need the galaxy centers, sky background values, 
and the orientation of the galaxy relative to the sky plane, estimated from the shape of the galaxy’s 
outer isophotes. In principle, this additional input for decompositions are published for all S^G 


galaxies in Munoz-Mateos et al. 2014, available via IRSA/S^G Pipeline 3. However, at the time 
our Pipeline 4 decompositions were made, these data was not yet available. Therefore, we made 
our own sky background estimates and ellipse fits. Also, the automatically created masks (see 


Munoz-Mateos et al. 2014) were visually inspected and hand-edited when needed (these edited 


masks later became part of the final S^G masks). If the decompositions are rerun starting from 
the output files provided in P4, it is important to use the data and mask images, as well as the 
pre-defined parameters offered via P4. 


In summary. Pipeline 4 consists of scripts for editing the masks, determining the galaxy centers, 
estimating the sky background, fitting isophotal ellipses, preparing the input files, and running 
GALFIT. It also includes tools for visualization of the GALFIT output files (see Sect |2.4[ ) and 
routines for storing the data on IRS A server (Appendix A) and the P4 web pages (Appendix B). 


2.2.2. Mask images 


The raw masks for the S^G 3.6 /im images were made in P2 with the SExtractor software 


(Bertin & Arnouts 1996), as described in Munoz-Mateos et al. (2014). Various automatic 

















detection thresholds for point sources were used. However, it soon became evident that no single 
criterion was sufficient to exclude all extra sources, without sometimes affecting also the galaxy light 
itself, in which case the masks needed manual editing. Also, in some cases the images contained 
artifacts that needed to be removed by hand. To speed-up this editing process, we developed a 
small portable IDL-routine (MASK_EDIT). Basically it displays on the screen simultaneously the 
original and masked images, and allows the user to remove/insert masked regions interactively. As 
an initial step of P4, all of the raw 3.6 /im masks were visually checked and edited if needed. The 
resulting masks are suitable for the purposes of our structural decompositions. However, because 
the wings of the PSF are quite extended (see Section 2.2.6) more extensive masking might be 
required in some applications. The MASK_EDIT routine, with source code and examples of use, is 
available at the P4 web page. 


2.2.3. Galaxy centers, sky background, isophotal profiles 


After the edited masks were completed, we run the galaxies through a semi-automatic IDL 
script which determines the galaxy centers, sky background levels and galaxy orientation param¬ 
eters. The accurate galaxy center is measured with the cntrd-routin^ after its approximate 
location is interactively defined. We also have the option to mark the center by force, in case the 
automatic center finding routine does not work satisfactorily even after repeated trials. 


The regions used for estimation of the sky background are identified manually, by selecting 
several (typically 10-20) locations outside the visible galaxy, while avoiding the image edges or 
contaminated areas. The local sky values in these locations are obtained by taking medians of the 
non-masked pixels in 40 pix x 40 pix boxes. The global sky background (SKY) and its uncertainty 
(DSKY) are then estimated from the mean and standard deviation of these local values, respectively 
(see Fig. [^. In section using the estimated DSKY, we show that the expected uncertainty of 
decomposition parameters caused by possible uncertainties in background subtraction is negligible. 
We also determine the average RMS sky variation, by taking the median of standard deviations in 
different sky regions (after removing outliers by iterative 3-sigma clipping). We use the sky RMS 
estimates in Section 2.2.5| for assessing the validity of theoretically calculated sigma-images. 


We calculate the isophotal profiles with a pyraf script called from IDL, using the standard 
IRAF ellipse algorithm ( Jedrzejewski [1987 ). As inputs for the ellipse fitting the sky background 
subtracted data image and the edited mask image are used. We fix the ellipse center to the 
previously found galaxy center and use a logarithmic increment of 0.02 between isophote levels. As 
often happens with IRAF ellipse, the fit does not necessarily converge over the whole galaxy area: 
we have an option to re-try the fit with different starting locations until a successful fit is obtained 


^cntrd is part of IDL Astronomy Library (Landsman 


1993). It locates the position where the brightness gradient 


IS zero. 
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over the whole galaxy region (see Fig. [^. From the isophotal profiles, we choose a semi-major axis 
range from which outer orientations {{h/a)outer = 1 ~ Pouter•> PPouter) are estimated. Also, a rough 
estimate of the galaxy outer radius, RgaU is made to define the image region used in the GALFIT 
decomposition. 


Figures and give examples of typical plots produced during these preparatory steps, il¬ 
lustrating the sky background fitting and the elliptical isophote profiles. The estimated eouter and 
PAouter are marked. In all our decompositions, we fix the orientation of the disk component to 
these outer value^ and interpret them to represent the galaxy viewing inclination. Therefore, 
extra care is taken to estimate the orientations reliably. For example, the corresponding incli¬ 
nation idisk — a)outer Is visually checked by de-projecting the galaxy images to face-on. 

Figure]^ shows an example of such a de-projection, also comparing the estimated inclinations with 
those calculated from axial ratios given in the HyperLeda database. Typically, our {hjcPjouter and 
PAouter are determined at much lower surface brightness levels than those in HyperLeda (which 
are mainly from RC3 de Vaucouleurs et al. ( 1991D ). P4 values are thus less affected by bulges, bars, 
or prominent spirals, and should reflect better the orientation of the underlying extended disk, 
which appears with circular outer isophotes in face-on projectioi][^ In some cases, S^G images are 
so deep that the outermost isophotes are dominated by an outer stellar halo rather than the disk. 
Good examples are NGC 681, NGC 1055 and NGC 4594. Possible misinterpretations of the outer 
isophotes were avoided by visually inspecting all the images: when the disk (identified with spiral 
arms, rings, and lenses) is clearly more inclined than suggested by the outermost isophotes of the 
image, the isophotes in the disk region were used for the estimate of galaxy orientation. For nearly 
face-on galaxies, the possible stellar halos are more difficult to distinguish, but in these cases the 
involved error in the orientation is less important. The final P4 axial ratios and position angles, 
center locations, and sky background values are listed in Table I. For each galaxy we also include 
a flag indicating the inclination uncertainty: ’ok’ indicates that outer isophote axial ratio should 
give a reliable estimate of idisk^ ’u’ indicates that the inclination is uncertain, while ’z’ indicates 
that the galaxy is close to edge-on (the axial ratio is not used for an inclination estimate). 

A scatter plot of P4 axial ratios versus HyperLeda values is presented in Fig. (upper left 
frame; only galaxies with flag=’ok’ are shown). As anticipated, the P4 axial ratios are on the 
average closer to unity than those in HyperLeda, though the difference is not very large (median 
( 6 /a)p 4 — {b/a)HyperLeda = 0.024). On the other hand the standard deviation of the difference is 
quite large (^ 0.1). The upper right frame makes a similar comparison to P3 axial ratios (Muhoz- 


Mateos et al. 2014) which correspond to a fixed surface brightness level /is .6 = 25.5 mag/arcsec^ 


On average, P4 orientations are measured at about 0.9 times this distance. The scatter is now 


^The reason is to reduce the degeneracy of different model components in decompositions 

^This expectation is of course not valid for a vertically extended (say T < 0) galaxy disk, nor in the case 
of intrinsically non-circular disks. However, the htted GALFIT expdisk-function assumes an inhnitesimally thin 
intrinsically axisymmetric disk, so any other treatment would be inconsistent in the decompositions. 
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significantly rednced and no systematic difference is seen between P3 and P4. The lower frames in 
Fig. I compares the position angles (now also ’z’ galaxies are included). In general the differences 
between P4 and HyperLeda are fairly small: for bja < 0.8 the median absolute difference is 2°. 
The difference between P4 and P3 is even smaller (median absolute difference 0.9° for bja < 0.8). 
Nevertheless there are some exceptions, most notably NGC4594 (The Sombrero Galaxy): in this 
case the P3 fixed isophote orientation corresponds to the extended halo, while the P4 orientation 
refers to the edge-on disk. 

Since we are fixing the disk orientations in the decompositions it is important to check the 
consistency of our inclinations. Figure displays the histogram of the P4 axial ratios for Hubble 
types —3 < T < 10. In case of a randomly oriented sample of thin disks, the distribution of b/a 
should be fiat. In case of finite vertical thickness a drop would be expected near a lower limit 
b/a = where qi is the intrinsic aspect ratio of the galaxies. According to Fig. [§i such a drop 
is evident for b/a < 0.15. However, overall the sample contains an excess number of galaxies with 
small axial ratios b/a < 0.5 (see the dashed line in Fig. [^). Similar trend is seen also when using 
the HyperLeda axial ratios (Fig. [^) or P3 isophotal orientations (Fig. [^). A possible explanation 
for the excess of small b/a ratios is that the S^G sample has been selected (Sheth et al. 2010) 
using an inclination-corrected blue magnitude limit {BTcorr = 15.5): if this dust correction were 
exaggerated, say for very late types, it would lead to an excess of faint, highly-inclined galaxies. 
This explanation is supported by the solid curves in Fig. which display the histograms when 
limiting to galaxies with non-corrected BT <15: now the histogram of P4 values is quite fiat. The 
histogram for P3 isophotal axial ratios is rather similar, though there are somewhat fewer small 
5/a < 0.2 values. This could be due to the above-mentioned faint stellar halos: in case of nearly 
edge-on galaxies a fixed surface brightness level could pick up the rounder faint outer envelopes, 
whereas in P4 we have in such cases tried to trace the disk isophotes. On the other hand, compared 
to both P3 and P4, the HyperLeda distribution has a clear deficit of large axial ratios, mostly likely 
due to the infiuence of inner non-axisymmetric structures. 

It is interesting to compare our sky background estimates to those in P3. In P3 an automatic 
sky measurement is made using 45 sky regions with 1000 pixels each. The regions are chosen close 
to the distance 2i?25 from the galaxy center (i ?25 is the blue band 25 mag isophotal radius from 
HyperLeda; if needed the distance of sky regions is modified manually). According to Fig. there 
is a very good agreement in the estimated sky backgrounds between P3 and P4 (see the right frame 
which takes into account that different PI mosaics are used for some of the galaxies). This good 
agreement is remarkable as the measurements are made completely independently and with different 
methods. The median difference between the sky determinations (0.0006 MJy/sr) is only about 1% 
of the typical sky background value, and its standard deviation (0.003 MJy/sr or ) is comparable to 
the magnitude of global sky variations in both sets of estimates (see Fig. [^. However, Fig. [^also 
reveals some cases where the difference between P4 and P3 is significant: inspection of the images 
indicates that this is due to a bright star (NGC1055), a nearby interacting component (NGC3327, 
NGC4647), or a too small FOV (NGC2655). In two cases (NGC1300, UGC10288) the final PI 


(Sheth et al. 2010 
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mosaic used by P3 is much improved over the earlier version used in P4. 


Fig. 0 compares our sky background variation estimates (’DSKY’ denotes global variations 
between sky measurement regions and ’RMS’ the average of the locally determined rms-scatter) 
with the corresponding estimates in P3 ( Munoz-Mateos et al. |[^014 ; their parameters ESKYl and 
SSKYl, respectively). There is a good overall agreement in the level of estimated global variation 
(left frame): the somewhat larger values for P4 are likely to follow from the larger range of radii 
we used for the sky measurement regions compared to P3. Also the local sky rms values show good 
agreement (right frame). 


2.2.4- Input data images 


As input for the GALFIT decompositions we use the 3.6 /im images. Because all necessary data 
reduction and calibration were already done in PI, the main preparatory steps are to subtract the 
estimated sky background value and determine which image region to include in the decomposition. 
In principle, GALFIT can also fit the sky background. However, this requires that the decomposed 
image region contains sufficiently large regions free of galaxy light or other contaminants. Use of 
such large image regions would slow down the decompositions considerably. Even more importantly, 
the S^G images often fill a substantial part of the raw frames or there are sudden jumps in the 
background levels (well outside the galaxy). To have a control of where the sky level is estimated. 


we chose to do the sky background evaluation manually, as described in Sect. |2.2.3| and to limit 
the decomposition to the rectangular region around the galaxy center. In practice, we choose 

Rfit = 1.3 X Rgai, where Rgai is our visually estimated outer size of the galax}|^ Finally, the image 
header keyword EXPTIME is set to 1 sec (as a default GALFIT will normalize the input data 
values with EXPTIME, which keyword is not relevant for PI mosaics), and all NaN^s (bad image 
values indicated with Not-a-Number value) are replaced with a constant value, and flagged in the 
mask in order to prevent them from affecting the decompositions. 


2.2.5. Sigma-images 

The sigma-images quantify the statistical uncertainty of each image pixel and thereby deter¬ 
mine the weights applied in GALFIT decompositions. This uncertainty contains two contributions: 
the noise contribution associated with the number of photons arriving at the instrument (’photon 
noise’ or ’shot noise’), and the noise originating from the instrument itself. The photon noise is 
assumed to follow a Poisson distribution, and it arises from two sources, the flux associated with 


®Later comparison to P3 isophotal radii published in IRSA indicates that the median < Rfn region/R 25.5 >= 1.7, 
where R 25.5 is the Pipeline 3 isophotal radius at = 25.5. The region is thus large enough to ensure that 

also the fainter outer parts of the galaxy are included in the ht. 
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the galaxy light and the flux coming from the sky background (zodiacal light). The main concern 
in the construction of the sigma-images is that the relative contributions of the photon noise and 
the instrumental noise are correctly estimated, so that correct relative weights are used in the 
decomposition for the bright central regions of the galaxies and for their faint outskirts. 

The sigma-images are calculated using the pixel values and header information in the 3.6 /im 
data images and the pixel values of the weight images. The images provided by PI are in flux units 
(MJy/sr), and for the calculation of the noise their pixel values F are converted to the number of 
electrons 


N,= 


F + Fi)a 


^ Fframe ^ frames ^ 9^ 


( 8 ) 


where F^g is the zodiacal light background which has been subtracted from the frame prior PI by 
the automatic Spitzer pipeline (its value is given by the header keyword SKYDRKZB). Note that the 
flux F contains besides the galaxy light also the sky background which has been subtracted in P4, 
F = Fgai + Fsky. The Fconv is the conversion factor between flux units and original digital units 
(header keyword FLUXCONV, in units of MJy/sr per DN/sec), T frame is the integration time/frame 
in seconds, Nframes is the number of combined frames for each pixel, and g is the detector gain 
factor (GAIN in units of e/DN). The number of frames combined is coded to the pixel values W 
of the weight images, Nframes — kP/10. Note that Tframe — 30 sec must be used instead of the 
original integration time/frame given by the header keyword FRAMTIME: this is because during the 
compilation of PI mosaics the pixel values have been normalized to this value regardless of the 
original The statistical uncertainty of N^, in each pixel is then calculated as a combination 

of Poisson noise (photon noise) and the readout noise of the detector (RON), 

(T^(iVe) = iVe + Nframes X R0N‘^. (9) 


We use RON = 15.0, 14.6, and 21 electrons, for FRAMTIME = 12, 30, and 100 secs, respectively. 
Note that these values, communicated by the Spitzer Science Center Helpdesk, deviate slightly from 
those given by the image header keyword RONOISE. The cr{Ne) is then converted to the estimated 
uncertainty of the image flux (note that a{Fgai) equals cr(F) since Fgj^y is constant) 

(yest{F) = cr(iVe) ^ Fconv! {Fframe ^ frames ^ ^)* (1^) 


In order to assess the validity of this estimate we compare it to the actual noise measured 
directly from the image. In Fig. this is done for the sky measurement regions. In the left 
frame the measured sky RMS (an average over all sky determination boxes) is plotted against 
the estimated acst from Eq. Colors distinguish between archival images from the cryogenic 


^This concerns the treatment of archival images observed during the cryogenic mission phase; all warm mission 
S^G observations have Tframe = SOsecs. 
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mission phase (original exposure time/frame either 12, 30, or 100 secs) and the new observations 
during the warm Spitzer mission (time/frame 30 secs, with the total exposure time of 240 seconds). 
For the archive images the overall agreement is quite good: there is a practically linear trend 
RMS ^ 0.9 (Jest holding for all three frame times, with the largest noise levels corresponding to the 
shortest frame times which have the largest contribution from the readout noise. The factor ^ 0.9 
is probably due to the PI mosaicking process, during which the images have been combined and 
sampled to 0.75 "pixel size from the native pixel size of 1.2 ". Because of this sampling the adjacent 
pixel values are strongly correlated, which is not taken into account in our theoretical estimate. 
Instead of trying to account in detail for the noise propagation during the mosaicking process we 
apply an empirical correction 


(^decomp{F) = 0.9 aestiF) (cryogenic mission) 


( 11 ) 


to be used in decompositions of cryogenic phase archival images. 

In contrast, for the warm mission the observed RMS is nearly 50% larger than the theoretical 
estimate (Fig. |^, indicating the presence of an additional source of noise. Also, there is a noticeable 
drop in RMS/(Jest ratio near the ecliptic plane (not present in the data from the cryogenic phase), 
indicating that the photon noise contribution to the (Jest (largest at / 0°) is overestimated 

compared to the instrumental contribution (constant with /). Following the advice of Spitzer 
Science Center Helpdesk, we include an additional instrumental noise component {(Jeonf)^ which is 
added quadratically to the theoretical noise estimate. To account for the PI mosaicking process, 
the multiplicative factor of 0.9 is again included. We thus adopt 

^decompi^) = 0.9 (j‘^est (F) + ^conf (warm mission) (12) 

for the warm mission images. The value of the empirical correction term (Jeonf estimated by 
this formula when applied to the sky measurement regions; the same formula is then applied to all 
image pixels. The adopted values of (Jeonf sire listed in Table I above. 

Figure illustrates the magnitudes of different contributions to the sky background noise. 
For the archival images (cryogenic phase, left frame) the noise is dominated by the readout-noise, 
though the Poisson contribution due to zodiacal light (the Fgpy we have subtracted in P4 + the 
SKYDRKZB subtracted during automatic Spitzer pipeline) still has a noticeable contribution. 
For the warm Spitzer mission (right frame) the extra noise term is even larger than the readout 
contribution. 


Nevertheless, at the central parts of the galaxies the photon noise associated with the galaxy 
light Fgai is the largest source of noise. This is illustrated in Fig. 10) which compares the Poisson 
and background contributions as a function of surface brightness. The two horizontal lines indicate 
the typical background noise levels for the cryogenic (lower) and warm (upper) phases (includes 
both instrumental and noise due zodiacal light). The inset figure illustrates how the a-map looks 
for the galaxy NGC3992 (observed during the warm mission). Near the center (/i 17), the 
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photon noise due to Fgai completely dominates, though already in the bar region (/x = 20 — 21) 
both photon and instrumental contributions are important. Altogether the sigma-images and thus 
the applied relative weights between galaxy and background regions are intermediate between those 
typically encountered when decomposing ground-based optical and NIR-images. In the former case 
the photon noise due galaxy light usually dominates, while in the latter case the sigma-image is 
almost completely dominated by the background noise, so that the weight is almost constant for 
all pixels (this applies e.g. to |Janz et ah 2014 GALFIT decompositions of Virgo dEs based on 
ground-based H-band images). 

In principle, the obtained (Jdecomp is just a statistical estimate of the true underlying variance at 
each pixel. We did some experimentation by smoothing the sigma -images (median averaging with 
kernels amounting up to 20 pixels). Except in the case of a few galaxies with very centrally peaked 
light profiles, this smoothing had very little influence on the final decomposition parameters. For 
the galaxies where smoothing played a role, the derived parameters were in any case uncertain (for 
example, the bulge Sersic index obtained unrealistically high values > 10). In the end, we decided 


to apply no smoothing at all. Tests related to the sigma-images are presented in Section 4.3 


2.2.6. PSF-image 


The IRAC data is not very well sampled: its native pixel resolution is 1.2", which is close to 


the Gaussian spread of a point source observed at channel 1. As discussed in detail in |Peng et al. 
I (2010), in such a case an oversampled PSF should be used. The IRAC PSF has also wide wings 


(see Fig. 11), so that a relatively large convolution box size must be used in decompositions: we set 
this to 40" X 40" (in some cases with a very centrally peaked light profile this region was extended 
to 150" X 150" with considerable increase in CPU time). Note also that IRAC PSF depends slightly 
on the instrument orientation. Therefore, in principle a separate PSF should be used with each 
image, determined from point sources in the same frame, or a combination of appropriate PSFs, in 
case the final image is a combination of several images obtained at different times. Clearly, such a 
procedure would be very time consuming. Fortunately, such an accuracy is hardly needed in our 
decompositions. The common oversampled PSF provided by T. Jarrett was used for all images, 
made as a composite over several instrument rotation angles. Fig. [redisplays the PSF, as well as 
a Gaussian profile (with FWHM — 2.1") approximately matching the core of the composite PSF. 
Also shown is an azimuthally averaged profile of the composite PSF. It will be shown in Section 
2.5 that it is important to account for the central core, as well as for the nearly circular wings 
of the PSF, whereas the outermost spikes have less importance for the obtained decomposition 
parameters. 











Table 1. 


Pipeline 4 parameters: galaxy center, outer orientation, and sky background 


IDE 

Xc 

Yc 

PA ± dPA 

ELL ± dELL 

RANGE 

FLAG 

SKY 

DSKY 

RMS 

O'conf 

ESOOll-005 

780.84 

299.02 

42.7 ±0.1 

0.747 ±0.002 

13 - 16 

z 

0.0125 

0.0039 

0.0109 

0.0092 

ES0012-010 

760.68 

294.95 

146.2 ±1.1 

0.542 ±0.012 

75 - 90 

ok 

-0.0022 

0.0025 

0.0109 

0.0096 

ES0012-014 

782.57 

459.99 

31.0 ±5.2 

0.580 ±0.040 

67 - 75 

u 

0.0080 

0.0033 

0.0106 

0.0091 

ES0013-016 

478.49 

282.78 

-14.3 ±3.1 

0.343 ±0.031 

52 - 82 

ok 

0.0041 

0.0027 

0.0101 

0.0084 

ES0015-001 

291.43 

294.80 

125.7 ±2.7 

0.586 ±0.018 

56 - 75 

u 

0.0050 

0.0018 

0.0102 

0.0086 

ES0026-001 

795.50 

495.47 

19.3 ±21.1 

0.060 ±0.028 

52 - 60 

ok 

0.0125 

0.0024 

0.0105 

0.0091 

ES0027-001 

776.34 

294.63 

12.5 ±10.1 

0.216 ±0.061 

112 - 127 

u 

0.0109 

0.0032 

0.0112 

0.0097 

UGC12791 

295.88 

289.63 

82.6 ±1.4 

0.709 ±0.037 

45 - 60 

ok 

0.0406 

0.0033 

0.0110 

0.0094 

UGG12843 

290.20 

287.90 

17.5 ±3.9 

0.553 ±0.058 

52 - 67 

u 

0.0418 

0.0037 

0.0105 

0.0086 

UGG12846 

571.15 

851.20 

-4.7 ±14.0 

0.127 ±0.031 

60 - 67 

u 

0.0446 

0.0007 

0.0019 

0.0000 

UGG12856 

291.07 

296.89 

16.9 ±1.5 

0.615 ±0.049 

48 - 75 

u 

0.0401 

0.0025 

0.0111 

0.0099 

UGG12857 

284.91 

280.12 

33.5 ±0.2 

0.760 ±0.009 

15 - 30 

z 

0.0453 

0.0018 

0.0110 

0.0096 

UGG12893 

294.93 

283.15 

87.2 ±3.8 

0.131 ±0.039 

52 - 67 

ok 

0.0463 

0.0022 

0.0110 

0.0093 


Note. — Galaxy center Xc,yc is given in pixels, ELL zb dELL and PA EdPA are the outer isophote ellipicity and position angle 
together with their standard deviations in the measurement range, given by RANGE (in arcsecs). EL AG indicates whether the 
inclination can be reliably estimated from the ellipticity {idisk = cos“^(l — Couter))- ok = reliable, u = uncertain, 2 ;= nearly edge-on 
galaxy. SKY^ DSKY, and RMS give the estimated sky level and its global and local variation (in MJy/sr). The last column (Jeon/ 
gives the estimated extra instrumental noise component during the Spitzer warm mission (See Sect 2.2.5). 



SKY= 0.04506 +/- 0.00211 


NGC3992.phot.1 .fits (2011 -03-31T21:04:07) 
NGC3992.1 .finmask.fits 


XC,YC(iraf)= 1316.68,920.88 
ELLIP = 0.428 +/- 0.025 
PA = 70.1 +!- 3.4 
range= 270.- 370. 

SKY= 0.04506 +!- 0.00211 

A25j3ix=322.9 T= 4.0 B 
BT= 10.46 MABS=-21.24 
ELLIP, PA= 0.302 77.5 

17-Aug-2011 10:23:27.00 



make_datalist0.pro (version HS230610/211210) 


NGC3992 



500 


-500 


-500 0 500 




Fig. 1. — Example of the determination of sky background for NGC 3992. The small upper left 
frame illustrates the center location found with cntrd-routine (black cross) relative to the image 
isophotes near the center (dashed red lines indicate the nearest integer pixels). In the lower left 
frame, red boxes indicate the local regions used for estimating the sky background: the mean and 
the RMS of these local median values were adopted for the sky background and its uncertainty 
(SKY and DSKY, respectively). The dashed ellipse indicates the visually estimated galaxy size 
{Rgal)- The white specks indicate masked stars. The upper right frame shows the intensity 
profile after subtracting the SKY value (indicated in the title of the plot; note the linear scale, 
intensities are from IRAF ellipse fits), marking also the ^DSKY (vertical red/green lines) and 
the median sky values in local measurement regions (boxes). The vertical dashed line corresponds 
to Rgai The insert shows the same profile, but as magnitude versus where a is the isophotal 

major-axis distance: a de Vaucouleurs profile would appear a straight line in this plot. The lower 
right frame shows the intensity profile in magnitude units (AB-magnitudes): red/green profiles 
correspond to adding/subtracting DSKY to the sky background. All distances are in pixels (0.75"). 
Similar plots for all sample galaxies are given on the P4 web site. 

The Figure + caption modfied according to referee’s suggestion 
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<PA>= 70.1 +/- 3.4 



NGC3992.phot.1.fits 



-400 -200 0 200 400 



Fig. 2.— Example of the isophotal profiles derived for NGC 3992 using the IRAF ellipse routine. 
The plots in the left display the PA and ellipticity profiles versus semimajor axis (a) of isophote 
ellipse: the dashed vertical lines indicate the range used in estimating the outer disk orientation 
parameters; solid horizontal line indicates the mean over that range. The upper right plot shows 
the observed, sky subtracted image, clipped at 1.3 x Rg^i (the image region used in decompositions). 
In the lower right the isophotes are plotted on top of observed (masked) image: the blue isophotes 
correspond to the a range from which the outer disk orientation {couter^ PPouter) was derived; the 
blue line indicates the assigned PAouter- Similar plots for all sample galaxies are available on the 
P4 web site. 
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NGC3992 



-200 0 200 


HYPERLEDA: pa,inc= 77.5 45.7 



-200 0 200 


P4: pa,inc= 70.1 55.1 



-200 0 200 


Fig. 3.— Example of deprojections with HyperLeda and P4 orientation parameters {{b/a)outer and 
PAouter) for NGC 3992. The red and green ellipses on the original image in the left illustrate the 
HyperLeda and P4 parameters, respectively (semimajor axis of the ellipse equals Rgai)^ the frame 
in the middle shows the deprojection with HyperLeda orientation, and that in the right with P4 
parameters: here the inclination is taken simply as i = cos~^(h/a) (note that this differs from the 
inclination (“incl”) listed in HyperLeda, which includes a morphological type-dependent correction 
for the assumed disk thickness). Clearly, the face-on disk is closer to axisymmetric when using 
P4 parameters: the difference would remain if the thickness-corrected HyperLeda inclination were 
used (47 degrees for this example). Similar plots for all sample galaxies can be found in the P4 
web site. 
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0.0 0.2 0.4 0.6 0.8 1.0 

b/a (P4) 



Fig. 4.— Comparison of outer disk orientation parameters to those in HyperLeda, and to the 
P3 orientation parameters corresponding to fixed /is.6 = 25.5 mag/arcsec^ isophote. In the upper 
row axial ratios are compared: the green line indicates the running median of HyperLeda (or P3) 
axial ratio, calculated in bins of 100 galaxies; gray indicates the RMS scatter in the bin. Red line 
indicates unit slope. The labels give the median difference and rms of the difference compared 
to P4: P4 axial ratios are generally larger than those in HyperLeda while the difference to P3 
is small. The lower frames displays the absolute difference in the position angles: squares mark 
deviant points with \APA\ > 15° for bja < 0.5. In the upper frames only galaxies with orientation 
uncertainty flag ’ok’ are included, while in the lower frames also galaxies with flag ’z’ are included. 
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0.0 0.2 0.4 0.6 0.8 1.0 

axial ratio (b/a) 


/home/heikki/P4_JUTTU_PRO/P4FIGjDainc_compareJeda_histo 


heikki@pc091079 Thu Sep 18 18:06:05 2014 


Fig. 5.— Distribution of axial ratios, calculated with P4 (a), HyperLeda (b), and P3 orientation 
parameters. Hubble types T < —4 and T > 10 are excluded (using Buta et al. 2014 Mid- 
IR classification). Dashed line indicates the whole S^G sample, with the magnitude selection 
BTcorr < 15.5, where BTcorr is the inclination-corrected blue magnitude from HyperLeda. Solid 
line corresponds to a similar limit, but using non-corrected blue magnitude BT. To ease the 
comparison the P4 histogram for BT < 15 is shown as the shaded region in each frame. 
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Fig. 6.— The left frame displays the used P4 sky background estimates (SKY) in comparison to P3 
(Muhoz-Mateos et al. llMii their variable ’SKYl’). The large apparent differences are due to use 
of different versions of PI mosaics. In the right we have corrected the P4 values to corresponds to 
the final PI mosaics (=those used in P3). Excluding the few deviant cases (discussed in the text) 
the median difference between P3 and P4 is 0.003 MJy/sr. Light and dark blue symbols indicate 
observations during cryogenic and warm Spitzer missions, respectively. 
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Global sky variation Local sky variation 



Fig. 7.— Comparison of P4 and P3 sky variation estimates. The left frame displays the global 
variations, estimated from the standard deviation of sky measurements at different areas (variable 
’DSKY’ in P4 and ’ESKYP in P3; the symbol colors are the same as in the previous Figure). In 
the right the local variation, estimated from the median scatter of sky values in local measurement 
areas (variable ’RMS’ in P4 and ’SSKYl’ in P3). 
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Fig. 8.— Comparison of measured and theoretical a. The left frame displays the measured 
RMS in sky regions against the estimated a, which takes into account the read-out noise, and the 
Poisson noise due to sky background (including the zodiacal light contribution removed by the 
automatic Spitzer pipeline). For the cryogenic mission phase the agreement is fairly good, with 
RMS ^ 0.9 (Jest (indicated by the dashed line; solid line indicates a one-to-one correspondence). 
For the warm phase the observed RMS is about 50% larger than the theoretically estimated noise. 
The right frame shows the ratio of the observed and estimated noise as a function of ecliptic latitude. 




0.03 




Fig. 9.— Different contributions to the sky background noise. The symbols indicate the measured 
average noise in sky measurement regions (RMS) in the cryogenic mission archival images (left 
frame; for clarity only those with 30 sec original frame time are shown) and in warm mission images 
(right frame), as a function of ecliptic latitude. The lines indicate various noise contributions and 
the total noise calculated with Eqs. (11) and (12): for clarity a mean over 20 deg bins is shown. 
Note that in the left frame the peaking of readout noise contribution close to ecliptic plane is just 
a spurious effect. 
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1 ^ 3.6 (AB) 


Fig. 10.— Comparison of galaxy flux and background noise contributions. The solid curve is the 
calculated Poisson noise associated with the galaxy flux Fgau as a function of surface brightness (in 
3.6 micron AB magnitudes; Fgai is converted to surface brightness with Eq. 13). The horizontal lines 
indicate the typical background noise levels for the warm and cryogenic mission phases (dot-dashed 
and dotted lines, respectively; they include both the noise associated with sky background flux and 
instrumental contributions). The insert shows the si^ma-map for NGC 3992. The structure in the 
background is due to different number of frames covering each pixel. Also notice how the galaxy 
stands out clearly on the sigma-map. The dotted line crossing the horizontal lines at /is.e ~ 25.5 
indicates the galaxy flux in MJy/sr. 
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Fig. 11.— The composite PSF used in decompositions (shown in the left frame in log-scale to 
emphasize the wings: in the right frame the PSF is normalized to its maximum value). The red 
solid curve in the right indicates a Gaussian with the same FWHM=2.1", approximately matching 
the inner part of the actual PSF (black symbols). The white curve on top of black symbols indicates 
an azimuthally symmetrized version of the composite PSF. The dashed blue curve indicates the 
cumulative flux outside a given distance from the center: the solid and dashed circles in the left, 
with radii 3.7" and 8.3" indicate the regions where 90% and 98% of the flux is concentrated (50% 
of flux falls within 1.3"). The composite PSF, oversampled with a factor of 5 (pix-size 0.15 arcsec, 
total size 30" x 30") was provided by Tom Jarrett. 
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2.3. Generation of input files for GALFIT decompositions 


The (ascii) input file for GALFIT specifies the galaxy data, mask, sigma, and PSF fits-files, 
and the region of the data image used in the decomposition. It also lists the components/functions 
used in the decomposition model, the initial guesses for the parameters, and specifies which of the 
parameters will be kept fixed, and which are iteratively varied in order to minimize the After 
convergence to a final solution, the final parameter values are written into an output file, with 
similar format as the input file. If needed, this output file can thus be used as an input for a new 
iteration (see Peng et al. 2002, 2010 for details). 


The input file also specifies how to convert the image values to magnitudes. The data images 
from PI are in flux units (MJy/sr). A conversion from pixel values Fi to (AB) surface brightness 
and integrated magnitudes is done with the formulas: 


/i3.6 = -2.5logiQ Fi + blogiQpix + zp 

magss = -2.5 logio 

i 


(13) 

(14) 


where pix = 0.75" and the zeropoint at 3.6 pm is zp = 21.097 (P3, [Mnnoz-Mateos et al. 2014). 
Values of pix and zp are inserted into GALFIT input file. 


All P4 input files for 1-component (Sersic) and 2-component bulge-hdisk (Sersic-hexponential) 
decompositions were generated automatically. Similarly, template files were created for the multi- 
component decompositions, which contained, in addition to bulge and disk components, entries for 
a Ferrers-bar, and a central unresolved PSF-component. The user then manually choose which 
components are fit and which functions used in the final model (see Sect for more details). In 
all our decompositions we keep the centers of the components fixed to the galaxy center. The 
cases were this is clearly not appropriate (galaxies with off-center bulges and bars) are noted in the 
parameter files. 


1. In 1-component input files initial guesses are needed for five free parameters: the Sersic 
index n, the effective radius i?e, the total magnitude m, the isophotal minor-to-major axial ratio 
g, and the position angle PA. The starting values of m and were taken directly from the data 
(total galaxy magnitude and half-light radius), for the Sersic index n = 2 is inserted as an initial 
guess, and q and PA were set to arbitrary values (0.9 and 10°, respectively). We thus avoided 
using the measured outer isophotes, to force GALFIT to search through a wider parameter space 
while minimizing the Typically 1-component fits converged after 10-20 iterations. When the 
fit did not converge, or if the final parameters were nonphysical (say, n > 10, g < 0.05, very large 
or small i?e), a new decomposition was started manually with new initial guesses. Usually this did 
not lead to any improvement, indicating that GALFIT is indeed very efficient in avoiding spurious 
local minima. 


2. The 2-component bulge-disk models apply a Sersic-function for the bulge: they thus need 
guesses for the same Sersic parameters as before, except that now these refer to the central compo- 
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nent. Accordingly, we used the initial guess (bulge) = 0.5 xR^ (image) and mi)uige = '^zma^e + 1- 
For the disk we use either ^expdisk’ or ’edgedisk-function^ depending on the estimated galaxy in¬ 
clination. In case of low or moderate inclination b/a ^ 0.2 (corresponding to i < 80°), we use 
the ^expdisk’ function, which needs two free parameters, the scale length {hr) and the integrated 
magnitude of the disk, mdisk- We chose hr — 0.25 xRg^i and rridisk — '^image + thus starting 
with a model with fairly massive and extended bulge. The disk orientation was fixed to the shape 
of the outer isophotes determined from the image (see Sect 2.2.3). In case of a nearly edge-on disk, 
b/a < 0.2, we use the ^edgedisk’ function with four free parameters: the central surface brightness 
/io, radial scalelength hr^ vertical scalelength /i^, and the position angle of the disk. The first guesses 
are ixo{disk) = /xo{image) -\- 3, hr as for the expdisk-model, while hzjhr = 0.1. The position angle 
is left free, with PAouter an initial guess. 


3. In the template files for the multi-component fits the initial bulge and disk parameters are set 
as for the 2-component models. For the Ferrers-bar the free parameters are the surface brightness 
at the effective radius of the bar, /ig, its outer truncation radius (denoted with Tout in Eq* 

(7)), its axial ratio, and its position angle. As initial guesses we choose /ig (bar)=/ig (image)-h3, 
Rbar = 0.25 X Rgah Qhar = 0.5, and PAi^ar = P^disk + 90°. For the magnitude of the unresolved 
central component we used m^sf = 'bUimage + 5. However, in practice we typically modified these 
pre-inserted template values even before starting the search of the final model, for example by 
adopting the output parameters from 2-component decompositions for the disk and bulge. 


2.4. Visualization of GALFIT decompositions: GALFIDL Package 


In its standard use, GALFIT is executed from the operating system command line, with 
an input file argument. This input file lists the input data files and the initial guesses for the 
parameters, as described above. The final decomposition parameters are written to an output file 
with a fixed name galfit.NN, where NN is a running number. Optionally, GALFIT makes a fits 
file containing the clipped data image (OBS; includes the region chosen for the fit), and total PSF- 
convolved model (MODEL), and the OBS-MODEL residual. Another GALFIT option is to write 
a FITS file containing model components in separate fits extensions. 


In P4 we have used GALFIT via GALFIDL, which is a set of IDL routines designed for easy 
visualization of the output from GALFIT decompositions. In addition, GALFIDL includes wrapper 
routines for calling GALFIT from inside IDL, with the advantage that the GALFIT output files 
and the produced plots are automatically renamed in a systematic fashion, using the names of 
the input files. We have utilized this by coding the galaxy identification and decomposition model 
components to the name of each produced output file (see Appendix A) 


The visualization options in GALFIDL follow those of the BDbar-decomposition program we 
developed earlier for the NIRSOS survey (Laurikainen et al. 2005), the most central of which is 
displaying a 2D plot of surface brightness vs. distance from the galaxy center (see Fig. 12). The 
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advantage of this, compared to the more commonly used azimuthally averaged profile, is that the 
contributions of different model components, with different apparent ellipticities, are easily high¬ 
lighted (Laurikainen et al. 2005; see also Gadotti 2008). The other visualization options include 
OBS-MODEL residual plots, profile cuts along a constant PA, comparison to observed profiles 
along isophotal major axis produced by IRAF ellipse, and plots showing schematically the different 
components included in the decomposition. The next section illustrates our decomposition strate¬ 
gies in more detail, concentrating on 2D-profiles. Additional plot types are illustrated in Appendix 
B, which describes the output released through a P4 web page for all S^G galaxies. 
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Fig. 12.— Schematic example of the 2D-profile plots. The black dots in the upper left frame show 
a synthetic observational image (exponential disk with added Poisson noise): the surface brightness 
(in arbitrary units) is plotted as a function of distance along the sky plane. The green dots indicate 
the best fitting GALFIT expdisk model; the width of the wedge-shaped profile depends on the disk 
inclination (here i — 60°). Note that this plot might leave an impression that the outer disk is 
not properly fitted: this illusion is due to the noise distribution appearing skewed in magnitude 
plots. The middle left frame illustrates the same model, after realistic noise (measured from the 
synthetic image) has been added also to the model profile, giving a visual confirmation that the 
model indeed is successful. The images in the right indicate the fitted model without (upper frame) 
and with noise (middle). The lowermost frames display a fit to a synthetic galaxy (image in the 
right) composed of a Sersic-bulge (here n = 2, with circular isophotes), an exponential disk, and a 
Ferrers bar. 
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3. Building the final multi-component decompositions - examples 


The final decompositions for S^G galaxies were done by fitting a maximum of 4 components. 
Typically the components were the bulge (5), disk (denoted either as D or Z, depending on whether 
the galaxy was close to edge-on), bar {bar) and the nucleus (A^, but could be any combination of 
these. The ingredients of the model are indicated by concatenating the designations of the compo¬ 
nents to the final model name: this same naming convention is used in the names of decomposition 
output files stored to IRS A (Appendix A). Note that the component designation is based on the 
intepretation of the component, not the function used in the fit. A disk (’D’), though most often fit 
with the expdisk function (1969 cases), may also be fitted with ferrer2 (69 cases) or sersic functions 
(99 cases). Similarly, in six cases a bulge (’B’) was fitted with an expdisk or edgedisk^ and in one 
case a ’bar’ with a sersic function. All elliptical galaxies were fitted with a single Sersic and are 
designated as B). 


In all final decompositions the orientation parameters of the outer disk were fixed and we also 
fixed a and /3 in the Ferrers function (<^=2, /3=0). All other parameters were left free for fitting. 
However, to find the structure components properly it was convenient to temporarily fix many of 
the model parameters at the beginning, and then release them one by one. For some galaxies, the 
length of the bar was kept fixed even in the final model. This was the case if GALFIT persistently 
gave a clearly incorrect bar length when compared to visual evaluation (in such a case the 
minimization was attempting to fit some other feature than a bar). 


Altogether over 20 different combinations of components were used in the final decomposition mod¬ 
els; Table collects an inventory of the main categories. This diversity of models is motivated by 
our desire to measure the bulge (if present) and the underlying disk parameters in a reliable man¬ 
ner. Note that our definition of ’bulge’ is quite broad, based on the excess flux in the central parts 
of the galaxy over that associated with the disk and bar components (’photometric bulge’). The 
decompositions themselves do thus not attempt to judge the physical character of this component, 
whether a merger-related, velocity-dispersion supported classical bulge, or a rotationally-supported 
’pseudo-bulge’ (Kormendy 1982), representing either a secularly formed central stellar disk com¬ 
ponent (Kormendy 1993) or a bar-related inner boxy/peanut component formed via bar vertical 
buckling ( Combes fc Sanders| 1981; [Athanassonla [2005[ ) . However, in Paper 2 we address the 
deduced bulge parameters (n, B/T) in the context of often-used classical/pseudo bulge indica¬ 
tors ( Kormendy fc Kennicntt||2004[ ) and also make comparisons to compilations of pseudo-bulges 
identified based on their HST morphology and star formation properties ( Fisher fc Drory||2010 ). 


In (non edge-on) galaxies with two distinct disk components (desgnated with dd, the inner disk 
was fit either with an exponential or a Sersic function, depending on the flattening of the profile. 
Such inner disk components differ from our photometric ’bulges’ by their much shallower profiles; 
they are also usually associated with a distinct inner spiral structure. Small central components 
were fit with the PSF, indicated as W’ in the model names. However, because of the limited 
resolution of S^G images, many of those structures, particularly in late-type spirals, might actually 
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be small bulges rather than nuclear point sources. Typically these components contribute less than 
a few percent of the total flux. 


3.1. Non-barred galaxies: 


The decompositions were made starting from simple 1 and 2-component models, and then 
adding as many components as necessary. For non-barred galaxies the process leading to the final 
model was the following: 


(1) Accepting the automatic 1-component model (single Sersic) as the final model. This was 
the case for elliptical galaxies (see NGC 3962 in Fig. [T^ . 

(2) Accepting the automatic 2-component bulge/disk decomposition as a final model. A typical 
example is NGC 3938 (Fig. [l4| ). 

(3) Adopting a bulge/disk model, after interactively finding modified initial parameters that 
converged to an acceptable final fit. 


(4) Adding a nucleus component or an inner disk (e.g. NGC 1357, Fig. 15) to the bulge/disk 
model. 


(5) When the galaxy had no obvious bulge we started from a single exponential disk, and if 
necessary, a second disk and/or nucleus was added (see NGC 723, Fig. 16). 


When the outer profile was affected by a possible stellar halo, the outermost part of the profile 
was not fitted. The best model was vetted by looking at the original image, the residual image after 
subtracting the model, the 2D surface brightness profile, and the ellipticities of the structures. The 
value of final was not used as a criterion in assessing the relative merits of the models (often 
a simpler final model was preferred even if a more complicated model would have yielded slightly 
smaller reduced 


Also, it is well known that many elliptical galaxies have small inner disks (Rest et al. 2001), 


and it has been shown by Huang et al. (2013) that many elliptical galaxies are better fitted with 


multiple Sersic profiles. Nevertheless, such a detailed approach was not taken in this study, in which 
the emphasis is in the analysis of disk galaxies (paper 2). It is worth noticing that while using deep 
images like those in in an automatic fit the bulge profile even in late-type spirals is easily 

degenerate with the outer part of the disk. In automatic fits this may lead to an unrealistically 
large Sersic n and for the bulge, of which NGC 1357 is a good example (Fig. 15). 
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Table 2. Main categories of final decomposition models. 


Disk: moderate inclination 



1889 


BD 

3II 



BDbar 

213 



ND 

214 



NDbar 

184 



Dbar 

458 



DD 

125 



D 

367 


Disk: nearly edge-on 



362 


BZ 

55 



NZ 

62 



Zbar 

8 



ZZ 

II3 



Z 

126 


elliptical: 

B 


26 

ALL 



2277 


Note. — Final decompositions were made for 2277 galaxies: in case 
of low or moderate inclination (apparent e < 0.8), the disk compo¬ 
nent was fitted with expdisk-iunction, while for nearly edge-on galax¬ 
ies ( e > 0.8), the edgedisk-function was used. In models BD and 
BZ, a bulge component was identified besides a disk, and it was mod¬ 
eled with a Sersic-function. These models may also contain additional 
disk components or unresolved central components (modeled with psf). 
The models BDbar include those bulge+disk systems which contained 
also a bar (modeled with ferrer2). In ND or NZ models the central 
component is modeled with PSF instead of ^Sersic-function. This may 
represent either a true central point source or (more commonly) an 
unresolved bulge. The models NDbar include also a bar. The models 
Dbar and Zbar have no inner Sersic or ps/-components, but include a 
bar-component. They may also contain an outer disk component. The 
DD models contain both an inner and outer disk (and no bulge nor 
bar), while D models refer to pure disks. Similarly Z models apply a 
single edgedisk-function., while ZZ models contain both thin and thick 
disk-components. 
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3.2. Barred galaxies: 


For barred galaxies a similar step-wise approach was followed. NGC 936 (Fig. © and NGC 
5101 (Fig. [Ts] ) are good examples demonstrating the importance of preventing the bar from mixing 
with the bulge flux. Adding a bar component to a simple bulge/disk model drastically changes 
the obtained properties of the bulge (for NGC 936 B jT drops from 0.46 to 0.19; for NGC 5101 
from 0.99 to 0.22). NGC 5101 has also a type II profile in the disk break/truncation classification 
associated to a broad outer ring (Laine et al. 2014). Using the edge of the ring as a manifesta¬ 
tion of a different flux distribution in the outer disk might be a bit misleading. Because of such 
ambiguities in the interpretation, we typically fit the type II disk profiles with a single exponential 
component. However, there are other barred galaxies in our sample, such as IC 4901 (Fig. 19), 
in which two exponential components (-h Ferrers function for the bar) were used for fitting the 
disk. In this particular galaxy using two exponentials is necessary, and those clearly correspond to 
distinct surface brightness components. Note that the outer disk of NGC 5101 is clearly lopsided 
(see the residual plot in Fig. 18). Such asymmetries are not taken into account in the pipeline 


decompositions (in case of strongly distorted galaxy no final model was made). See Zaritsky et al. 
(2013) for a detailed study of galaxy lopsidedness using S^G images. 


3.3. Pure disk galaxies: 


A third main group of galaxies in our sample are those having no obvious bulge. They may 


have a single exponential disk (NGC 3377A in Fig. 20), or more than one disk components (NGC 


723 in Fig. 16). The structure fit as an inner disk in NGC 723 consists of broad, prominent, and 


tightly wound spiral arms. Bulgeless galaxies may have bars, of which NGC 3517 (Fig. 21) is an 
example. To get an homogeneous estimate for the scale length of the disk for these galaxies, the 
outer disks were always fit with an exponential function, even in galaxies where the disk would 
have been somewhat better fit by a Sersic function with n slightly less than unity. Generally, the 
assumption of an exponential disk is good, but there are also cases, like ES0026-001 (Fig. [2^ in 
which a Sersic function would actually be a better choice. 


3.4. Edge-on galaxies: 


The GALFIT models for the nearly edge-on galaxies assume that the disk is viewed completely 
edge-on. A bulge, and in some cases also a bar or an additional thick disk component were included 
(Fig. [ 2 ^ . In these models also the vertical thickness was an output parameter. However, these 
models are tentative, and are meant solely as starting points for better, scientifically oriented 
decompositions. There already exists detailed modeling of edge-on galaxies in S^G, based on 


fitting their vertical profiles to hydrodynamical thin-thick disk models (Comeron et al. 2011 


2012). Their radial luminosity profiles have been analyzed in( Comeron et al. 2012; Comeron et 


















25 


al. 2014; Martm-Navarro et al. 2012). 


3.5. Scope of pipeline decompositions 


The P4 models for the spiral galaxies T > 0 are generally good, giving reliable estimates 
for parameters such as the bulge-to-total flux ratio (S/T), the scale length of the disk (hr) and 
its central surface brightness (/io)- However, despite the fact that up to four components were 
fit, the pipeline decompositions for the early-type disk systems (T<1), because of their complex 
structures, are often insufficient. These systems may have nuclear bars, ovals and lenses, which 
are not included in our models in any systematic fashion. Because of this, the pipeline B/T 
flux-ratios, particularly for SO galaxies, can be over-estimated. Including all these structures will 


require even more complex decompositions, such as those done in the near-IR by (Laurikainen et 


al. 2005, 2006, 2009, 2010). Such time consuming modeling goes beyond the scope of our current 


P4 decompositions; nevertheless, the P4 decomposition output files provide good starting point for 
further fine-tuning. 
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Fig. 13.— Figs. - [^ show the examples of final decompositions described in the text. In 
the big panel the galaxy image is shown in an inverted logarithmic scale (magnitude range 27 > 
/is.6 (AS) > 18), clipped to display the main morphological characteristics. The three small panels 
show the masked original image (upper right panel), the model image (lower right), and the residual 
OBS-MODEL in the middle (range ±lma^; white indicates excess light over the model). The lower 
left frame shows the 2D profiles of the observed and model images (black and white dots), together 
with the model components (colors; labels indicate the relative fraction of flux in this component; 
in this particular case there is only one component). The same components are also marked, with 
the same colors, on the lower right model image: the semimajor axis of the ellipse corresponds 
to 2i?e of the component in question. The mid-IR classification from Buta et al. (2014) is also 
indicated. In this particular example for NGC 3962, the single Sersic fit provides an acceptable final 
model. The overall profile is close to a de Vaucouleurs profile (Sersic n = 5.6) in accordance with 
the morphological classification (E). Nevertheless, the slight bends in the profile and the structure 
in the residual image suggests that if desired, it would have been possible to get an even slightly 
better fit by including multiple components. 








- 37 - 


NGC3938 SA(ls)c 
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Fig. 14.— NGC 3938: Example of a galaxy in which the automatic Sersic bulge + exponential disk 
fitting gave an acceptable final solution. In the profile plot the spirals appear as small undulations 
on the generally well-fitted disk. The labels in the profile plot indicate the relative fraction of flux 
in various components. 
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NGC 1357 (R’L)SA(LS)0/a 



Final BDD-model: n=4.1 




Automatic BD-model: n=10 
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Fig. 15.— NGC 1357: The surface brightness profile of this galaxy shows a small bulge and a large, 
fairly exponential disk. However, the automatic bulge-disk fit would give an unreliably large bulge 
extending through the whole galaxy (lowermost row). A more reasonable fit is obtained by adding 
another exponential disk component to the inner part of the galaxy (upper profile). This inner 
component corresponds to the region of tightly wound spiral arms with higher surface brightness. 
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NGC 723 (R’)SA(rs)bc 
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Fig. 16.— NGC 723: This is a disk galaxy with no centrally concentrated bulge component, but it 
is clearly not a single exponential galaxy. The bump in the surface brightness profile corresponds 
to the strong high surface brightness spiral structure in the inner parts of the galaxy, here fit with 
a nearly flat part (Sersic function with n = 0.12). 
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Automatic BD-model: B/T=0.46 
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Fig. 17.— NGC 936: An example of a barred galaxy where inclusion of a bar component to the 
model makes a large difference in the bulge parameters. The lowermost row shows the automatic 
bulge/disk model, whereas the upper panels include a bar component. In the simple model, the 
bar ffux is degenerate with the bulge ffux. Multi-component decomposition is thus essential for 
getting a realistic bulge-to-total ffux ratio: the B/T = 0.46 in the 511-model, but drops to 0.19 in 
the BDhar model . 
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NGC5101 (R,R’2)SB(rs)0/a 
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Automatic BD-model: B/T=0.99 
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Fig. 18.— NGC 5101: Another example of BDbar decompositions where inclusion of the bar 
components is essential in getting realistic bulge parameters. Note also that in the final model the 
underlying disk is fit with an exponential function, although the outermost profile is downbending 
(Type II break/truncation). However, the steeper outer slope seems to be associated with a broad 
double outer ring, rather than a fundamentally distinct outer disk component. 













- 42 - 


1C 4901 SB(r)bc 
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Fig. 19.— IC 4901: An example of a barred galaxy, in which the disk is fitted with two exponential 
functions. Additionally, a small central psf-component is included, marked as a cross on the model 
image. The inner disk corresponds to the higher surface brightness part of the disk outside the bar, 
where the spiral arms are also prominent. 
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NGC3377A Im/Sph 
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Fig. 20.— NGC 3377A: A bulgeless disk galaxy, well fit with a single exponential function. 
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IC3517 (R’)SAB(s)dm 
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Fig. 21.— IC3517: Another example of a bulgeless disk galaxy. The surface brightness profile 
is well fit with a single exponential function. However, the image shows also an elongated inner 
structure, which can be fitted with a Ferrers function. 
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ES0026-001 (R’ 2 )SAB(s)cd 
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Fig. 22.— ES0026-001: An example of a galaxy in which the disk is fitted with an exponential 
function, although a Sersic function with n <1 would have given a more precise fit to the disk. 
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Fig. 23.— ES0533-004: An example of an edge-on galaxy fitted with two disk components: the 
ratio of thick disk to thin disk scale heights is = 3.9 and the thick disk contains 35% of the 

total model luminosity. For comparison, the detailed vertical profile fits in Comeron et al. 2012 
indicated = 4.3 — 4.6 depending on the radial location, with about 45% of light in the thick 

dist component. 
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4. Uncertainties of the decomposition parameters 

The formal uncertainties of the decomposition parameters have little significance, as they 
refer to purely statistical uncertainty due to image noise based on the assumption that the model 
is accurately describing the true underlying light distribution. Taking into account the complex 
morphology of most galaxies, this assumption is clearly not valid (see Peng et al. 2010 for detailed 
discussion of errors)]^ Related to this, the final value of the reduced ^ poor indicator of 

the goodness of the fit (even for a good model it is typically much larger than unity) and is thus 
not used as a decisive factor in choosing the preferred final model. In practice, the choice of 
the final model components plays a crucial role: for example as seen in Section omission of 
the bar component when a bar is present may lead to seriously biased bulge parameters. In this 
Section we perform a systematic comparison of bulge and disk parameters between 2-component 
and final multi-component models. We also first examine the potential uncertainties related to the 
preparation of data before the decompositions, namely the used PSF-function, the effect of sky 
subtraction uncertainty and the sigma-image. 


4.1. PSF 


As illustrated in Fig. in the IRAC PSF has extended wings. Moreover, the PSF and the 
orientation of its asymmetric extensions vary from image-to-image, which has not been taken into 
account in our decompositions. To check the importance of the PSF wings, we compared differ¬ 
ences in decomposition parameters obtained when the adopted composite PSF was replaced with 
a Gaussian PSF having the same FWHM = 2.P'. Fig. 24 compares the resulting effect on the 
Sersic parameters in 1-component models. Clearly, decompositions with the Gaussian PSF yield 
n values that are systematically too small, differences reaching even tens of percents for some of 
the galaxies (though the median deviation is less than 5%). However, these rather large deviations 
are not representative of the true uncertainties, but rather give an idea of the magnitude of the 
error if the tails of the PSF were altogether ignored. A better measure of the actual uncertainty 
in P4 decompositions is obtained by comparing with an azimuthally symmetrized version of the 
composite PSF. Clearly, now the differences in n are much smaller (see the red symbols in Fig. 24). 


We also checked the influence that the PSF has on the multi-component models. For that 
purpose we rerun all final decompositions that included both bulge and disk components {+ possi¬ 
ble bar and center components; total of 524 models after excluding nearly edge-on galaxies), using 
both the Gaussian PSF and the symmetrized PSF. Tablelists the median of relative differences in 
bulge H/T, n, Rg? disk scale lengths /i^, and bar-to-total ratio Bar/T, when compared to the results 
obtained using the standard composite PSF. The largest differences are seen for the Sersic param¬ 
eters while using the Gaussian PSF, whereas hr is barely affected. On the hand, the differences in 


®The formal uncertainties calculated by GALFIT are listed in the headers of pipeline output hies in IRS A 
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parameters between those obtained using the composite PSF and using the symmetrized version 
are negligible. Based on these results we conclude that the spikes of the PSF have no significant 
effect as long as the nearly circular wings of the PSF are included. The use of single composite 
PSF for all S^G images should thus be acceptable. 


4.2. Sky subtraction 

In principle, poor sky subtraction can severely affect the decomposition results, in particular 
the parameters of the disk. To constrain the possible magnitude of such uncertainties, we rerun the 
multi-component decompositions that included both bulge and disk components (+ possible bar 
and center components; same 524 models as above). Two additional sets of sky values, SKY' — 
SKY ± DSKY^ were used, where DSKY was the standard deviation of the different sky regions. 
Fig. 1^ shows the effect on the scalelength of the disk. Although individual changes can in few 
cases be large (/i^(mod)//i^(ori) > 1.2 in 9 cases when too small a sky is subtracted), the median 
differences are less than 2% (and even smaller in the other parameters of interest, see Table |^. 
The sky subtraction is not a concern in the current decompositions. 


4.3. Sigma-image 


The weights applied to various pixels have an important role in decompositions, in particular 
when the galaxy structure is complicated, so that the differences between the applied model and the 
true structure are large. As mentioned in Section 2.2.5 the cr-image itself is a statistical estimate 
of the underlying a in each pixel, so it might be reasonable to smooth it before applying it in the 
decompositions. In Fig. (left column) we examine the effect of sigma-image smoothing on the 
derived bulge parameters. A median filter is applied with a width of 5 pixels. Clearly the effect is 
quite small except for a few deviant cases marked on the plot. In these cases the bulge parameters 
are sensitive also to changes in the PSF or the sky background level. 


For comparison. Fig (right column) also illustrates the changes in bulge parameters if a 
constant sigma is assumed at all image pixels. A constant sigma exaggerates the relative weight of 
the central regions compared to the outskirts. Besides a large scatter, also a systematic increase 
of the estimated n is obvious: the median nmodluori = 1.2 (the mean ratio is 1.4). What typically 
happens is that the fit tries to reproduce the central peak with an increased n, even if the outer 
disk then becomes too much bulge dominated. Indeed, the bias (and the scatter) is particularly 
large for earlier type disks (open circles in the plot indicate T < 5; median rij^odluori — 1-25). This 
comparison reminds us that when decomposition parameters from different studies are compared, 
it is also important to pay attention that similar weights have been applied. 
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4.4. Two-component versus multi-component decompositions? 

Automatic 2-component Sersic-exponential (or Sersic-Sersic) models are often applied to large 
data surveys. This is a natural approach as the data quality (depth/angular resolution) might be 
insufficient for more detailed modeling so that the large effort in multi-component decompositions 
does not seem justified. Moreover, it has been recently claimed (Tasca & White 2011) that 2- 


component decompositions (Sersic -h exponential) are sufficient also for fitting barred galaxies. 
Their argument was based on obtaining similar average BjT ratios for barred and non-barred 
galaxies in their 2-component bulge/disk models. They reasoned that if the omission of the bar 
were a problem it should manifest as a higher B/T for barred galaxies. However, to accurately 
address this matter one has to compare the different types of decompositions (2-component and 
multi-component) for well-defined samples of barred/non-barred galaxies. 


Such a comparison between different decomposition models is shown in Fig. ^7j Again, those 
galaxies for which the final model contains both a bulge and a disk are studied. For the non-barred 
galaxies (those with no 6ar-component; leftmost column in the Figure) the bulge parameters (Sersic 
n, H/T, Re/hr) in automatic 2-component runs are almost identical to those in the final models. 
This agreement is expected because over 80% of the final non-bar models are just Sersic-expdisk 
models (15% have two disk components, and 2% have an extra central component), and typically 
the automatically found BD models did not need any refinement. For barred galaxies (those with a 
6ar-component in the final model; middle column), the obtained median values depend drastically 
on whether the bar is included. This result emphasizes that the examples of decompositions given 
in Section]^ highlighting the importance of modeling the bar (e.g. Figs. [T7| and 18) were not 
exceptional cases. Overall, ignoring the bar increases the estimated B/T ratios by a factor of 
2-3 because of gross (even by as much as a factor of 5) overestimate of Re and n. For example, 
for spirals in the range 1 < T < 5 the 2-component decompositions suggest n > 4 whereas the 
multi-component runs indicate n ~ 1 — 2. Altogether, in the final models the difference in bulge 
parameters obtained in the multi-component decompositions for barred and non-barred galaxies is 


fairly small (right column in Fig. 27). 


The conclusion that multi-component decomposition models are essential to measure realistic 


2007; Gadotti 2008 


bulge parameters for barred galaxies is not new (Laurikainen et al. 2006 

2009[ )). A similar conclusion, based on synthetic images, was reached also by 


Weinzirl et al. 


Laurikainen et al. ( 2005D . 


In Fig. 28 we compare the combined bar/non-barred sample of the previous figure with the 
decompositions in Laurikainen et al. ( 2007| ). Because of the large fraction of barred galaxies, the 
difference in the obtained bulge properties between the 2-component and multi-component decom¬ 
positions remains significant, even when barred and non-barred galaxies are considered together. 
We find an excellent agreement between the current multi-component results and those in |Lau 


rikainen et al. (2007), obtained with a different decomposition code (BDBAR] however BDBAR 
uses IDL Curvefit and is thus based on the same Levenberg-Marquadrdt minimization as GALFIT), 











































- 50 - 


and based on different near-IR image data. It is worth noticing that in these decompositions the 
Sersic n for Hubble types Sa-Sc is nearly n ^1, whereas in the decompositions by |Tasca fc White 


(2011), for the same Hubble types, the Sersic index is peaked at n ^4. Small values of the Sersic 
index, similar to ours for these Hubble types, are reported also by Graham & Worley (2008). 
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Table 3. The effect of a modified PSF on final decomposition model parameters 


GAUSSIAN PSF 


SYMMETRIZED PSF 



median(D) 

median( D ) 

median(D) 

median( D ) 

BjT 

-1.5 % 

4.5 % 

-0.1 % 

0.2 % 

n 

-3.0 % 

7.8 % 

0.0 % 

0.5 % 

Re 

8.9 % 

9.8 % 

0.1 % 

0.3 % 


0.0 % 

0.3 % 

0.0 % 

0.0 % 

Bar/T 

-1.5 % 

4.4 % 

-0.1 % 

0.2 % 


Note. — D stands for the relative difference (e.g. D = [n(mod) — n(ori)]/n(ori)), where ’ori’ 
refers to the standard composite PSF. Medians are used to characterize the typical deviations and 
the scatter, to eliminate spurious cases where the decompositions with Gaussian PSF converged to 
a different type of solution. 
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Fig. 24.— The effect of the PSF on decomposition parameters. The plots on the left show the 
Sersic index and effective radius in 1-component fits using modified PSFs instead of the standard 
composite PSF: black symbols indicate results using a Gaussian PSF (wings truncated), and red 
points when using a symmetrized composite PSF. Larger black (red) symbols indicate points devi¬ 
ating by more than 25% (10%) from the unit line. In the right panels the histograms of the relative 
changes in the parameters are shown: black and red colors have the same meaning as in the left 
frames. 
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Table 4. The effect of sky subtraction on final decomposition model parameters parameters 



SKY + DSKY 
median(D) 

median( D ) 

SKY - DSKY 
median(D) 

median( D ) 

BjT 

0.2 % 

1.5 % 

-0.1 % 

1.6 % 

n 

-1.4 % 

1.8 % 

1.8 % 

2.0 % 

Re 

-1.2 % 

1.4 % 

1.7 % 

1.8 % 

h'p 

-2.5 % 

2.6 % 

2.6 % 

2.8 % 

Bar/T 

0.2 % 

1.5 % 

-0.1 % 

1.6 % 


Note. — D stands for the relative difference (e.g. D = {n(mod) — n(ori))/n(ori))^ 
where ’ori’ refers to the standard sky subtraction. 
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Fig. 25.— The effect of sky background subtraction on the disk scalelength. All hnal decomposition 
models including both a bulge and disk (and possibly additional bar and central components) were 
rerun using images where the assumed sky background was modified by ±DSKY, where DSKY 


is the conservative estimate of global sky variations in the image (see Section 2.2.3). 
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Fig. 26.— Sensitivity of estimated bulge parameters (Sersic n, bulge-to-total flux ratio S/T, 
effective radius Re in kpc’s) on the used a-image. In the left column, we have smoothed the P4 
(j-images with bpix x bpix median filter, while in the right it has been replaced with a constant a. 
The scatter plots show the modified parameter values versus the original ones. In the right, the 
red open and black filled circles refer to galaxies with mid-IR type T < 4 and T > 5, respectively. 
Lines corresponding to one-to-one correspondence are drawn in each frame: in the right uppermost 
frame the dashed line indicates rimodluori — 2. 
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Fig. 27.— Comparison of bulge parameters between automatic 2-component and final multi- 
component decompositions. The morphological type T is from the Buta et al. (2014) mid-IR 
classification. Comparison is made for the galaxies for which the final model included both bulge 
and disk components. In the left frames decompositions for non-barred galaxies are compared, 
while the middle frames show those with a bar component in decompositions. In the right frames 
the final decompositions for barred and non-barred galaxies are compared. The symbols stand for 
median values in bins with five or more galaxies, error bars are errors of the mean values in the bin. 
Note that the Buta et al. (2014) classification contains also half-integer values of T, resulting from 
averaging over two rounds of classification. However, the number of galaxies with half-integer values 
is much less than those with integer T. Therefore, when binning the galaxies we have rounded the 
half-integer values randomly to the nearest smaller or larger integer value; same is done in Fig. [33| 
below. 
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Fig. 28.— Comparison of final S^G decompositions parameters to those of Laurikainen et al. (2007; 
the values are tabulated in Table 2 in Laurikainen et al. 2010) multi-component decompositions for 
NIRSOS {Ks band, 143 galaxies with —3 < T < 1 ) and OSUBSGS (i/band, 129 galaxies 2 < T < 7) 
data. The S^G results show 524 galaxies for which the final decomposition model contained both a 
disk and a bulge component (excluding nearly edge-one galaxies). For comparison, also the results 
of semi-automatic 2-component decompositions are shown. The symbols stand for the median 
values in each bin with 5 or more galaxies, while the error bars denote the error of the mean. Note 
that here we use the optical classifications from HyperLeda, to facilitate comparison with previously 
published results. 
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4.5. Disk breaks 


One of the main goals of Pipeline 4 is to obtain measurements for the galaxy size-magnitude 


scaling relations. In order to be consistent with earlier analysis (e.g. Courteau et al. (2007)) the P4 
final models as a default use single exponentials for the disk. However, deep optical and near-IR 


surveys (Erwin, Beckman & Pohlen 2005; Pohlen & Trujillo 2006; Gutierrez et al. 2011; Muhoz- 


Mateos et al.|2013 ) have shown that only a fraction of galactic disks 1/3) are simple exponentials 
(=Type I in [Pohlen fc Trnjillo||2006 classification). Instead the typical brightness profiles consist of 
two (sometimes three) exponential subsections with different radial slopes. When the outer disk has 
a steeper slope, the galaxy is classified as possessing a Type II break (’truncation’), and conversely 
if the outer slope is more shallow, it is classified as a Type III break (’antitruncation’). [Kim et al 
(2014) have recently made 2D decompositions for 144 barred S^G galaxies taking into account disk 


breaks in their decompositions with the BUDDA code ( de Souza et al.||2004 ; Gadotti||2008 ; Gadotti 
| |2009| ). Their fitting function for the disk consists of two exponential sections, with different scale- 
lengths {hin and hout) inside and outside the break radius Rhreak- They also made decompositions 
where they fitted the disk with a single exponential component. Their result indicate that the 
inner scale lengths for two-component disks are typically about 40% longer than the scalelengths 
obtained in single disk fits; they thus conclude that “it is important to model breaks in Type II 
galaxies to derive proper disk scale lengths.” 

Nevertheless, it is not always obvious what is the ’’proper” disk scale length estimate to use in 
various scaling relations, in case the galaxy exhibits several exponential subsections. For example. 


it is well known (Pohlen & Trujillo 2006) that Type II breaks are often connected to outer rings 


associated with bar OLR resonances. Such breaks are indeed dominant for early type barred disks 
(T < 3; Laine et"ar]|2014 ). Since the bar torques are able to push material from the CR regions out 
toward OLR, this will promote a shallower distribution inside the break radius. However, beyond 
the OLR the effect of bar is insignificant, so that the underlying disk can remain more or less 
intact. In such a case it might in fact be the outer, rather than the inner scalelength that would 
better characterize the original overall mass distribution. On the other hand, for later Hubble 


types the Type H break is often connected with the end of prominent spirals (Laine et al. 2014) 


and could be due to suppressed star formation: for such a case the inner scale length might indeed 


be more appropriate to characterize the disk as a whole. Laine et al. (2014) also find that for such 
spiral-related breaks the ratio hinner/houter is typically closer to unity than for OLR related breaks. 


Figure compares Pipeline 4 decompositions with several recent disk truncation studies, 
which use subsamples of the same S^G dataset. In this plot the disk scalelengths are displayed 
against the stellar mass derived in P3 ( Mnnoz-Mateos et al. [[2014 ). Besides the above-mentioned 


Kim et al. (2014) 2D decomposition study, we also compare with Muhoz-Mateos et al. (2013) and 


Laine et al.j ( |2014 ), where fits to 1-dimensional profiles were conducted. First of all, the Figure 
(upper row) indicates a very good agreement for the scale lengths of Type I profiles between all 
four studies, conducted with independent methods. Secondly, it illustrates the significant difference 
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between the inner and outer slopes for Type II (and III) profiles, amounting to roughly a factor of 
two (see Mnnoz-Mateos et al.||2013 , Kim et aL~||2014[ ) . The P4 single disk scalelengths seem to fall 
quite close to being a geometric mean of hinner and houter derived in earlier studies. 


To emphasize the possible ’unperturbed’ nature of Type II outer disks, we compare in Fig. 
the P4 scale lengths vs. stellar mass for Type I galaxies with the type II outer scale lengths derived 
in the above mentioned disk break studies. Indeed the differences between the Type I single disk h 
and the Type II houter are quite small, much smaller than the differences compared to hinner- The 
fits to the data also give the impression that hinner I houter ratio gets closer to unity for less massive 
galaxies: this is in accordance with the above mentioned dominance of spiral-related less abrupt 
truncations for later and thus on average less massive spirals. 


The Pipeline 4 single exponential fits have a convenient feature of representing an effective 
average over inner and outer disks (when both present). They thus provide a homogeneous set 
of robust scale measurements, not sensitive to factors modifying the local slopes. Nevertheless, a 
possible caveat is that the fitted effective single h might become dominated by different degrees 
by the inner/outer parts, depending on the galaxy surface brightness. For example, the estimated 
h might be biased toward hinner when the disk central surface brightness decreases toward less 
massive galaxies: this would be the case if the image depth was not sufficient to cover the galaxy 
regions beyond the break radius. Fig. addresses this potential problem by comparing the trends 
of the break radii with respect to galaxy mass, to that of the galaxies’ visual outer extent {RgaU see 
Sect 2.2.3; a similar trend would result if i?25.5 were plotted instead of Rgai)> The figure indicates 
that a break, if present, should be detectable through the whole range of S^G galaxy masses. 


In summary, we feel confident that the single disk fits provide a useful overall estimate of 
the disk original scale length (and its extrapolated surface brightness), though especially in case 
of barred massive galaxies secular evolution might have led to significant deviations from simple 
exponentials, important to include in detailed models for individual galaxies. Moreover it is likely 
that the slope differences associated with breaks are smaller for later types, which form a vast 
majority of S^G galaxies. 


Nevertheless, as concluded by Kim et al. (2014), estimates of other decomposition parame¬ 


ters, such as the B/T ratio for massive galaxies would become more accurate if the inner slopes 
are accounted for (say, leading to less disk light assigned to bulge). The situation is somewhat 


analogous to the benefit of including additional inner components like bars (Laurikainen et al. 


2005; Gadotti||2008 ), lenses in SOs (Laurikainen et al. 2010), or barlenses ( Laurikainen et aL]|2014 ) 


into decompositions. However, for the goals of Pipeline decompositions, the expected magnitude 
of changes (about 10% relative change in B/T according to Kim et al.) is quite small, compared 
to the uncertainties related to choice of the decomposition model components (say, including a bar 
versus ignoring it). The choice of the code might also sometimes have a bigger effect. For example, 
Kim et al. ( 2014[ ) use NGC 936 as an example of Type II galaxy (see their Fig. 4). For this galaxy 
they fit a break at 98" and derive hinner — 53" and houter — 28", all very close to the measurements 
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in both Munoz-Mateos et al. | ( |2014 ) and Laine et al. (2014). On the other hand, the Pipeline 
4 single disk fit (see Fig. [TtI ) gives h — 40". We verified that truncating the disk in GALFIT 
decompositions at the break radius given in Kim et ah, reproduces their inner slope quite well (we 
get 57"). At the same time, the B/T we obtain increases slightly (from 0.19 to 0.22), as anticipated 
by Kim et alj^ Nevertheless, the B/T we obtain after accounting for the more shallow inner slope 
is still nearly 50% smaller than the value obtained by Kim et al. (S/T=0.32), probably because 
of some model/code dependent factors, such as how the image pixels are weighted, or the PSF is 
treated. 


^We also checked the effect of letting the boxiness and shape parameters of the bar free but these turn out to be 
small 
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Fig. 29.— Comparison of P4 single disk scalelengths with decomposition studies including disk 
truncations. In the left frames the hr vs. stellar mass, obtained in P4 (blue symbols) are compared 
with the results in Muhoz-Mateos et al. 2013| (M2013) and Laine et al. [2014 (L2014), where the 
inner and outer scalelengths (/imner and houter^ denoted with black and red symbols, respectively) 
were estimated for S^G galaxies from fits to 1-dimensional profiles. On the right, similar comparison 
to Kim et al. 2014 (K2014), who used 2D BUDDA decompositions for 144 barred S^G galaxies. 


In the uppermost frames Type I disks (no breaks) are compared: the lines show orthonormal fits 
to the measurements (orthogonal deviations minimized using the IDL PCOMP routine). In the 
middle same for Type II profiles: the dashed blue line indicates a fit to geometric means of the 
inner and outer scalelengths. In the lowermost frame: galaxies with Type III breaks. 
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Fig. 30.— Comparison of P4 scalelengths for Type I galaxies (blue points), with the outer scale- 
lengths of Type II profiles (red points; these combining Munoz-Mateos et al. 2013, Laine et al. 
(2014) and Kim et al. (2014) measurements). Also shown are the corresponding orthonormal fits 
(solid blue and dashed red lines, respectively). The black dashed line shows a fit to inner scale- 
lengths derived in the above mentioned studies (individual measurement points not shown). Note 
that the trend of single scalelengths in Type I galaxies resemble much more the outer scalelengths 
in Type IPs rather than the inner scalelengths. 
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Fig. 31.— Comparison of the break radii between the inner and outer disk segments in the studies 


of Muhoz-Mateos et al. (2013), Laine et al. (2014) and Kim et al. (2014); large and small circles 


refer to Type II and Type III breaks, respectively. Note that the break radius drops rapidly with 
galaxy stellar mass. Also shown by small squares is Rgah fhe P4 visual estimate of the galaxy 
extent in S^G images. The plot suggest that even in the case of low-mass, low surface brightness 
galaxies the depth of the S^G images is sufficient to assure that the break, if present, is not likely 
to be buried in the sky background. 
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5. Decomposition parameters 

In the current paper (paper 1) we provide all the l-component and final multi-component 
decomposition parameters in a tabular format (Tables 6 and 7), together with quality assignment 
flags. A brief check of how the major categories of the final models distribute among different 
Hubble types is also shown. All actual analysis will be presented in paper 2 . 


5.1. Quality assessment 


The full S^G sample contains 2352 galaxies, chosen according to their internal extinction 
corrected B-magnitude {Mscorr < 15.5), apparent B-band 25-mag isophotal diameter (D 25 > 
60"), galactic latitude {\b\ > 30°), and HI recession velocity {Vradio < 3000A:m/5), obtained from 
the HyperLeda database. Due to its mag-limited character, it contains a large number of low 
surface brightness late-type spirals and irregulars. Also, galaxies with peculiar morphology were 
not specifically excluded. In some cases the field-of-view (FOV) is not large compared to the galaxy 
size (the new Spitzer observation mapped regions covering at least 1.5 D 25 , but this condition was 
not fulfilled by all the archival galaxies included in the sample). In such cases the sky background 
is difficult to estimate reliably, and for galaxies near the ecliptic, the background may have larger 
gradients (see Muhoz-Mateos et al. 2014). Altogether, the sample contains a number of galaxies 
for which decompositions are less reliable, or not possible at all to carry out. 


Because it is important to estimate the reliability of the derived structural parameters, we 
have assigned to each galaxy a quality flag, running from 1 (worst case) to 5 (most reliable). The 
judgment was done partly by visual inspection of the original data and partly by evaluating the 
decomposition models: 


• Quality=l (31 cases) 

Reasons: Bad original data (very bright overlapping stars, strongly varying background, image 
defects). 

Action: Excluded from all analysis: galaxy identifications are listed in parameter tables but no pa¬ 
rameter values are given. P4 web page illustrates the raw data -h mask, but not any decomposition 
models. 

• Quality=2 (44 cases) 

Reasons: Original data is more or less fine, but the FOV is too small for reliable sky estimation. Al¬ 
ternatively, galaxies exhibit strongly distorted shapes which make even l-component fits unreliable 
(mergers, interacting, peculiar, strong warp, very lopsided). 

Action: l-component Sersic fit was done and the parameters are given in Table (and in the 
web pages), together with a comment indicating that they need to be taken with caution. Multi- 
component decomposition parameters are not given. 





-64- 


• Quality=3 (61 cases) 

Reason: Original data is fine but the galaxies have complex structures that require detailed multi- 
component models beyond the scope of the pipeline decompositions (which have a maximum of 4 
components). 

Action: 1-component Sersic fit is considered reliable (Table . Multi-component decomposition 
was also made and parameters are listed in Tablewith a cautionary comment. Web pages show 
both 1-component and multi-component decompositions. 

• Quality=4 (406 cases) 

Reason: Original data and decomposition are of good quality. However, the galaxy was either 
highly inclined (contained a ’z’ component; 333 cases), or it had complicated structure, so that 
there might be a degeneracy between model components (such as between inner and outer disk 
components) (73 cases). 

Action: All decomposition parameters given in Tables andand illustrated in the P4 web-pages. 
However, these are omitted from the analysis of disk central brightness and scale length in paper 

2 . 

• Quality=5 (1810 cases) 

Reason: Both the original data and the decompositions are of good quality. 

Action: All decomposition parameters are given in Tables and and illustrated in the P4 web¬ 
pages. 

Table summarizes the number of galaxies in different quality categories. 


5.2. One-component fits 


The output parameters of 1-component Sersic fits are listed in Table For 1-component fits 
the parameters are the Sersic index n, effective radius i?e, integrated magnitude mag^ axial ratio 
g, and major axis position angle PA (the centers are fixed to those given in Table I; the isophotes 
are assumed to be elliptical and to have fixed a shape and orientation with radius). Additionally, 
there is a column indicating the reliability of the fit. 


Single-component Sersic-fits are routinely used in large data surveys. This gives objective, 
easily reproducible results, that provide useful characterization of the galaxy global characteristics. 
For example, |Cappellari et al. ~| ( |2013 ) argued that Sersic n>4 largely finds the most massive early- 
type galaxies (M > 3 • IO^^Mq), which are also the slow rotators in their kinematic classification. 
Also, there are well-known correlations between galaxy color and Sersic index. Therefore, in paper 
2 we will present detailed analysis of the 1-component Sersic parameters for the S^G sample. Here 
we report just the dependence of n on the morphological type. 


Figure [32] displays the histogram of Sersic index-values in the 1-component models. Galaxies 
are divided into three bins according to their mid-IR Hubble type (E with T < —4, SO with 
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—3 < T < 0, spiral with 1 < T < 9, and irregulars T = 10). Clearly, the distribution of galaxies 
peaks at n 1 — 2, with a broad tail to larger values of n. There is also a secondary peak close 
to n ^ 4 (corresponds to a de Vaucouleurs profile), but this is not very prominent. The overall 
distribution reflects the nearly exponential profiles of many late type spirals, which dominate the 
S^G sample. For irregulars, the distribution peaks at n ^ 1. For SO’s, the distribution is much 
broader. The distributions remain essentially the same for less inclined galaxies, for example if the 
sample is limited to those with apparent bja > 0.5. 


5.3. Multi-component decompositions 


For multi-component decompositions the tabulation is more complicated, as model compo¬ 
nents/functions vary from one galaxy to another. Also, the same function can be used to de¬ 
scribe different structure components in different galaxies. We have decided to present the multi- 
component parameters in two different formats, one that is compact and easily human-readable, 
another more suited to automatic reading. In the first format (Table the first entry for each 
galaxy indicates the used model and the number of components. The next lines, for each component 
included in the model, indicate the physical interpretation of the component {B, D (or Z), bar, A), 
and the GALFIT function used (sersic, expdisk, edgedisk, ferrer2, psf), followed by the component 
parameters. The Table caption specifies which parameters are listed for each function. The other 
table (available via IRSA and P4 web page) lists for each galaxy all possible components and their 
parameters: empty values indicate that this component was not included in the decomposition of 
this galaxy. 

The S^G sample contains 358 galaxies which were considered to be close to an edge-on view 
and are excluded from further analysis in paper 2. Also, 26 are elliptical systems, modeled with a 
single Sersic function. This leaves 1855 moderately inclined disk systems. As discussed in Section 
there are over 20 different combinations of functions/components used in the final models, so that 
there is a need to group the decompositions in to major categories. A natural approach is to base 
this grouping on whether the decomposition model has a bulge component. Because the bulge can 
be modeled either with a “sersic” or “psf” component, depending on its apparent size, we have two 
categories, BD and ND models, respectively. When there is no trace of a bulge, the system can be 
either a single disk (D), possess a bar-like component (Dbar), or contain both inner and outer disk 
components {DD). 


The numbers and relative fractions of galaxies in these categories as a function of Hubble 
type are shown in Fig. 33 Here the mid-IR classification from Buta et al. (2014) is followed. 


Apparently the relative fraction of Ri^-models decreases gradually towards the late-type spirals. 
However, taking into account that most of the ND models describe small bulges (support for this 
claim is given in paper 2), indicates a much smoother distribution of galaxies with bulges, dropping 
rapidly only above T ^ 5. The fact that ND models cover a relatively large range of Hubble types, 
including SOs, is quite interesting because it indicates that SOs can possess very small bulges. This 
is in agreement with Laurikainen et al. (2010), where the same conclusion was made based on 
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decompositions of NIRSOS data. The result is consistent with the idea that at least some SOs might 
be former late-type spirals with small bulges, devoid of gas, following quenching of star formation. 


A very small BjT ratio in an early-type spiral has also been reported also by Kormendy & Barentine 


[ (2010). Another interesting feature in Fig is that many galaxies that lack bulges, can still have 
bars. These Dbar galaxies peak at Hubble types T — 7. Beyond T — 9 they are replaced with 
single disks, becoming almost the sole type of models for irregulars (T=10). The HH-models are 
most common (about 10%) for T = 9. 

Finally, Fig. gives examples of galaxies in these major categories, indicating the model 
components for 4 low mass and 4 large mass systems in each category. For BD and ND categories 
the barred and non-barred models are also distinguished. Note that ’barred’/’non-barred’ refers 
here to whether or not a bar-component was included to the the final decomposition model, not to 
any morphological classification; for example a non-barred BD model has been adopted for NGC 


5985, which has a SAB family classification (see Fig. 33). A detailed comparison to Buta et al. 
(2014) classification will be presented in paper 2. 


Table 5. Summary of decomposition quality flags 


Quality 

# 

#(>quality) 

1-component 

Multi-component 

Disk fiQ and hr 

1 

31 


- 

- 

- 

2 

44 

2321 

uncertain 

- 

- 

3 

61 

2277 

ok 

uncertain 

- 

4 

406 

2216 

ok 

ok 

uncertain (or z) 

5 

1810 

1810 

ok 

ok 

ok 


Note. — Quality flags (1-5) assess the reliability of decomposition parameters. The second 
column indicates the number of galaxies in each class, while the third column indicates the 
number of galaxies with decomposition parameters of this or better quality. 











Table 6. Parameters of 1-component Sersic fits 


Identification 

Gomment 

mag 

q 

PA 

n 

Re 

ESOOll-005 


14.52 

0.231 

43.46 

1.329 

18.61 

ES0012-010 


13.31 

0.493 

156.97 

1.592 

56.38 

ES0012-014 


14.67 

0.411 

23.92 

0.884 

47.35 

ES0013-016 


12.75 

0.508 

168.64 

1.314 

38.85 

ES0015-001 


14.39 

0.357 

110.22 

1.131 

34.34 

ES0026-001 


12.52 

0.689 

58.26 

3.126 

62.93 

ES0027-001 


10.75 

0.765 

74.38 

5.423 

96.07 

UGC12856 


14.04 

0.316 

18.48 

1.339 

41.22 

UGC12857 


12.99 

0.195 

33.50 

1.316 

19.09 

UGG12893 


13.58 

0.861 

94.74 

1.179 

33.55 


Note. — Comment indicates cases where no fit was made, fit failed to 
converge, or its parameters are not reliable, mag is the the total magnitude, 
q is the axial ratio and PA the position angle of elliptical isophotes, n is 
the Sersic index and Re the effective radius (in arcsec). 
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1 -component Sersic-fit 



Fig. 32.— Histogram of Sersic-index n in 1-component Sersic-fits. Dashed line indicates the 
whole S^G sample, after elimination of 61 galaxies due to bad image quality, bright nearby star 
etc. Additionally, for 8 galaxies the decomposition failed to converge (n = 20 in the plot). The 
main morphological types are shown separately with different colors: Ellipticals ( T < —3), SOs 
(—3 < T < 0), Spirals (0 < T < 10), and Irregulars (T = 10). Additionally there are 8 galaxies 
without classification (T = 99; the median value of their n — 2.3), and 34 dwarf galaxies (T = 11; 
median n = 1.3) which are not shown. The Hubble stage T is from the mid-IR morphological 
classification by Buta et al. (2014). The vertical dashed lines indicate the values n = 1 and n = 4, 
corresponding to exponential and de Vaucouleurs profiles, respectively. 





















Table 7. Parameters of final multicomponent decompositions 


#1 

ESOOII-005 

B 

Z 

_bz 

sersic 

edgedisk 

NGOMP=2 

0.571 

0.429 

quahty=4 

14.822 

20.912 

0.354 

43.902 

40.447 

9.764 

48.937 

2.955 

3.001 

#2 

ESOOI2-OIO 

D 

BAR 

_dbar 

expdisk 

ferrer2 

NGOMP=2 

0.944 

0.056 

quahty=5 

13.445 

22.177 

0.458 

0.390 

-33.795 

24.245 

33.953 

2I.I4I 

23.095 

#3 

ESOOI2-OI4 

D 

BAR 

_dbar 

expdisk 

ferrer2 

NGOMP=2 

0.835 

0.165 

quahty=4 

14.795 

24.704 

0.420 

0.249 

31.032 

II.102 

33.152 

83.009 

24.393 


#2352 

UGGI2893 

B 

D 

_bd 

sersic 

expdisk 

NGOMP=2 

0.021 

0.979 

quahty=5 

17.826 

13.617 

0.786 

0.869 

102.674 

87.208 

6.559 

20.150 

0.515 

22.134 


Note. — The first row for each galaxy is the running number (1-2352). The second row gives the galaxy 
name, the type of final decomposition model (coded to all output file names together with underscore- 
prefix), the number of components in the model, and the quality flag. If no final decomposition was 
made (quality= 1 or 2) then for this galaxy we set type=’-’ and NCOMP=0. The next NCOMP entries 
give: (1) the physical interpretation of the component (B-bulge, D-disk, Z-edge-on disk, BAR-bar, N- 
unresolved central component), (2) the GALFIT function used for it, and (3) the component’s relative 
fraction of the total model flux. The next entries depend on the GALFIT function. For sersic they are: 
mag^q^PA^n^Re^ for expdisk: mag^q^ PA^hr^ for edgedisk: go^PA^hr^hz: ior ferrer 2: go, q, PA, Ri^ar: 
and for psf. mag. Here mag is the total 3.6 gm AB magnitude, go is the central surface brightness in 
mag/arcsec^ (face-on brightness for expdisk and edgedisk, sky brightness for ferrers2), Re,hr^hz are in 
arcsecs. All decompositions assume a fixed common center for all components and elliptical isophotal 
shape, constant over radius. If there is an (outer) disk, its q and PA are kept fixed to those in Table I. 
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-5 0 5 10 

T 


Fig. 33.— Distribution of final decomposition model categories as a function of mid-IR Hubble 
stage from Buta et al (2014; the half-integer values of T have been rounded before binning, see 
caption of Fig. 27). In the upper panel the distribution of the original S^G sample (dotted 


histogram), and of galaxies with final models (red dashed line histogram; excluding the edge- 
on galaxies). The green histogram is the distribution of models where both ’bulge’ and ’disk’ 
components where identified (BD). This is also the subsample used in Section 4.4 when comparing 
automatic 2-component and final multi-component decompositions. The dark blue line shows the 
distributions for models with both ’nucleus’ and ’disk’ (excludes those with ’bulge’) (ND). The 
three other lines are for models with neither ’bulge’ nor ’nucleus’: ’Dbar’ stands for models where 
a Ferrers-bar was included (together with one or more ’expdisk’ components), ’DD’ stands for 
models with inner and outer disks, while ’D’ stands for a single ’expdisk’ model. In the lower 
panel the relative fractions of different models are plotted, normalized to the total number of non 
edge-on models (red dashed curve curve in the upper panel). Here ’BD’ and ’ND’ include both 
models with/without ’bar’: for clarity, these are not shown separately as the differences in the 
distributions are small. 
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Fig. 34.— a) Examples of different main types of final decomposition models. Four low mass and 
four large mass galaxies of each type are displayed. The left frames display the 3.6 /xm image, with 
fixed AB surface brightness range [18,27], while the right indicates the model components: the 
semi-major axis corresponds to 2i?eff of the component. The labels in the left frame indicate the 
galaxy name and the Buta et al. (2014) mid-IR classification; in the plots the Buta et al. underline 
notation is indicated with slanted characters. Labels in the right frames give the physical coding 
of the decomposition model components (same as used in the names of the decomposition files), 
and the logiQ{Mstar) (stellar masses are from Mnhoz-Mateos et al. |^014D . The colors of the 
ellipses indicate the used functions: expdisk (green), sersic (red), ferrer2 (blue). Similar plots for 
all galaxies are given in the P4 web-page. In a) examples of BD models (have bulge & disk but no 
bar component) are displayed. 
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Fig. 34D.— Examples of BDbar decomposition models (bulge &; disk with a bar component). 


’barf’ indicates that the length of the Ferrers bar was fixed in the decompositions. 
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Fig. [34|3 . — Examples of ND decomposition models (nucleus & disk, no bar). The central compo¬ 
nent (unresolved in decomposition) is indicated with a red dot. 
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34i.— Examples of NDbar decomposition models (nucleus & disk, with a bar). 
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343 .— Examples of Dbar decomposition models (disk & bar). 
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Examples of DD decomposition models (two disk components) 
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Fig. 


345 .— Examples of D decomposition models (single disk) 
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6. Summary and Conclusions 


Two-dimensional multi-component decompositions, using GALFIT3.0, have been performed 
at 3.6 /im wavelength for the complete S^G sample (2352 galaxies). Reliable decompositions were 
possible for 2277 galaxies. Quality flags are given for each galaxy based on our confidence on the 
model parameters. The main goal of the decompositions was to estimate the structural parameters 
of the bulges and disks in a reliable manner, which dictated our decomposition strategy. Most 
importantly, a bar-component was included in the decomposition model whenever present, to pre¬ 
vent its light from biasing the derived bulge and disk parameters. For the same reason, the models 
sometimes included a central point source and additional disk components. However, no attempt 
was made to match the detailed shape or length of the bar. 

We present automated single Sersic, 2-component bulge-disk (Sersic -h exponential) decompo¬ 
sitions, and human-supervised, individually checked multi-component models. In the final multi- 
component models, a maximum of four structural components are fit: bulge (Sersic), disk (exponen¬ 
tial), bar (modified Ferrers), and the nucleus (PSF). Different combinations of component functions 
were used. For example, in some barred galaxies it was convenient to fit the underlying disk with 
two different functions. As a first step, we estimated the sky background levels, derived the ori¬ 
entation parameters with ellipse fitting, and edited the masks to eliminate foreground stars and 
image defects. In general we found an excellent agreement with the independent P3 measurements 
(Muhoz-Mateos et al. lIMil ). 

The uncertainties related to the sky background, the adopted PSF-function, and to the treat¬ 
ment of sigma-images were tested. The decomposition data are released in IRSA, and in the P4 
web-page, where the decomposition models, the ellipse fitting, and sky background determinations 
are illustrated. The IDL-based tool (GALFIDL) used in visualization of GALFIT decompositions 
is also available on the web pages. Besides the decomposition output files, all input files needed 
in re-doing the decompositions are given in IRSA. All of this provides the possibility to refine the 
pipeline models for the needs of specific scientific goals. 


In particular, such refined models will be needed for early-type galaxies which often contain 
more structures than handled by current pipeline decompositions: for ellipticals such structures 
include nuclear point sources ( Lauer|[T985 ; Cote et al. 2006) and inner disk structures (Kormendy 


&; Bender 

1996, 

Buta et al. 

2014; modeled in Kormendy &; Barentine 

2010 

, Weinzirl et al.||2014) 


for SOs and early-type spirals the various lens structures ( Kor mend^| 1979 ; Laurikainen et al.|2009 ) 
should be accounted for, including the barlens components ( Laurikainen et aL]|2011 ) recently iden¬ 


tified as the more face-on counterparts of boxy/peanut bulges seen in nearly edge-on galaxies (Lau- 


rikainen et al.|[Ml4 ; Athanassoula et ^|2014| ). Such refined models become particularly important 


with the ongoing extension of S^G data to include 465 additional gas-poor early-type galaxies 
( Sheth et al.||2013 ), which were not part of the original sample which contained only galaxies with 
emission line velocity measurements. 


The main results are the following: 
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(1) Automatic single Sersic fits. The Sersic indexes peak at n ^1.5, having only a minor peak 
at n ^4, reflecting the fact that a large majority of the sample galaxies are spirals with extended 
disks. 

(2) Automatic 2-component bulge-disk decompositions: Such decompositions would suggest a 
large difference in the parameters of the bulges between barred and non-barred galaxies. Since this 
is an artifact caused by the inadequate decomposition model, we strongly caution against using 
simple bulge-disk decompositions for barred galaxies. 

(3) Final multi-component decompositions. In contrast to 2-component models, in our final 
models the differences in bulge parameters between barred and non-barred galaxies disappear, 
leading to the values of Sersic n ^ 1 — 2 for bulges in both types of galaxies. It means that if bars are 
not included in the fit, the flux of the bar is erroneously mixed with the bulge flux. This conclusion 
is consistent with several previous studies using a similar multi-component decomposition approach. 

(4) Small bulges containing at most a few percent of the galaxy flux at 3.6 fim appear in a 
large range of Hubble types, including SOs. This is in agreement with Laurikainen et al. (2010) 
where a similar result was obtained in near-IR. 

(5) At intermediate Hubble types (T=5-7) the very small bulges gradually disappear but the 
galaxies can still be either barred or non-barred. At the very end of the Hubble sequence pure disks 
become dominant. 

A detailed analysis of the properties of bulges, disks, and bars, as a function of morphological 
type, and other overall properties (galaxy mass and global color) will be presented in forthcoming 
paper (Paper 2 in preparation). 
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A. Appendix: Decomposition pipeline products in IRSA 


The results of pipeline 1-component and final multi-component decompositions are available 
via the IRSA database. For each galaxy, decomposition output parameters (outgal-Sle)^ and input 
fits-files are given. The user can refine/improve the given multi-component models by including 
more components or by utilizing additional GALFIT options for the component functions, 


10 


1) Input data for GALFIT decompositions 


IDE is the galaxy designation (e.g. NGC1097) 


fits-files: 

IDE.phot.l_nonan.fits 


IDE.1.finmask_nonan.fits 
IDE.phot.l_sigma.fits_ns 
PSF-1.composite.fits 


= 3.6 micron image used in decompositions. 

Bad pixel values (NaN’s) removed 
header modified to make GALFIT work correctly 
= corresponding mask-file 
= -"- sigma-image 
= PSF-image 


2) Output from GALFIT decompositions 


a) ascii-files: 


IDE_onecomp.outgal 
IDE_twocomp.outgal 
IDE_M0DEL.outgal 


Automatic best fit parameters for 1-component S\’ersic model 
-"- for 2-component sersic+expdisk (or sersic-edgedisk) fit 
Final decomposition model with up to 4 different components 


IDE is the galaxy designation (e.g. NGC1097) 

MODEL-string identifies the components included in the final multi-component model: 


b’ 

indicates 

bulge-component 

d’ 

indicates 

disk -"- 

bar’ 

indicates 

’barf’ -> 

non-axisymmetric structure, mainly bars 

length of the bar was fixed in the decompositions 

n’ 

indicates 

nucleus (or nonresolved bulge) 

z’ 

indicates 

edge-on disk 


is important to use the data and mask files from IRSA P4 directories when refining the given decompositions, 
instead of using corresponding data products from PI and P3 directories. Namely, the P4 outgal-^\es assume sky 
subtracted data images (with NaN image values removed) and EXPTIME keywords set to 1 sec. Additionally, the 
(P4 vs (PI & P3)) data and mask images may have small spatial shifts (a few pixels) and correspond to slightly 
different sky background levels, depending on when the various pipeline products were finalized. 
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e.g. "NGC1415_bdbar" -> bulge+disk+bad final decomposition model 
These outgal-files contain the decomposition output parameters 

Together with the input fits-files, the user can immediately repeat/refine the decompositions 
starting from the outgal-file (e.g. outgal -ol NGC1415_bdbar.outgal -> re-creates the decomposition) 


b) fits-file (for the FINAL model) 


IDE_M0DEL.outgal_subcomps.fits Final decomposition output images: 

extension 1 = OBS image 

extension 2,3,4... = model components 

File header contains also final decomposition parameters. 

c) jpg-files (for the FINAL MODEL) 

IDE_MODEL.outgal_profile.jpg Decomposition model compared with observations: 

- shows surface brightness at each image pixel vs distance from galaxy center 
Observed image, model image, and model components displayed separately 

- Collects also decomposition input & output parameters 
Labels indicate relative contribution of model components 


IDE_MODEL.outgal_residual.jpg Model-observation comparisons: 

upper row: clipped 3.6 micron image, masked image 
lower row: model image, OBS-MODEL residual 


IDE_MODEL.outgal_ldprof.jpg Decomposition model profiles compared with observations: 

- shows surface brightness as a function of isophotal semi-major axis, 
comparing IRAF ellipse fits to the observed image and to the 
model image (using isophotes of the observed image) 


IDE_MODEL.outgal_components.jpg Schematic plot of model components: 

- Upper row: Observed image and model image, with different model 
components marked: colors correspond to profile plots, and the semimajor-axis 
of the ellipse is 2 times the effective radius of the components 

- Lower row: same as upper row, except projected to the disk plane 
(assuming zero-thickness). Empty in case on edge-on final model. 


The final decomposition models, as well as various intermediate steps involved, are illustrated in the web-page 
\protect\vrule widthOpt\protect\href{http://www.oulu.fi/astronomy/S4G_PIPELINE4/MAIN}{http://www.oulu.fi/astronomy/S4G_PIPELI]) 
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B. Appendix: Decomposition Pipeline web-pages 


The input data used in decompositions, the various steps of the decomposition pipeline, and 


the final decomposition results are illustrated on the Pipeline 4 web sit6^^ 
http: //www. oulu. f i/astronomy/S4G_PIPELINE4/MAIN 
The web site consists of three layers of pages: 

1) Main page 

2) Index pages 

3) Decomposition pages 


The main page gives full instructions regarding the contents of the pages, followed by an alpha¬ 
betical list of all 2352 galaxies in the original sample Clicking on any of the galaxy names opens a 
corresponding index page^ with contains information for 100 galaxies, near and including the chosen 
galaxy (Fig. [3^ . 

The icons in the index page indicate the data available for each galaxy 

1) PI image mosaic 

2) clipped PI image, with mask 

3) clipped PI image 

4) Deprojected image 

5) Elliptical isophote profiles 

6 ) Galaxy center/sky background plots 

What decompositions available: 

7) 1-component Sersic-model (icon shows residual plot) 

8 ) Final-model (icon shows residual plot) 

9) Final-model (icon shows model components) 


An empty icon signifies no data/decomposition model, for instance if the galaxy was discarded 
from further analysis because of a nearby bright star, etc. In particular, the absence of a final model 
indicates that the galaxy was considered too problematic to fit reliably (e.g. closely interacting, 
very peculiar or warped). 


Clicking on the galaxy name links to its decomposition page (Fig. 36), which summarizes 
the mask, ellipse fitting, sky background, and galaxy center determinations, as described in Section 


The website does not contain actual fits-files, which are available via the IRS A server/P4. 
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In addition, the pages display for all three types of models (1-component Sersic model, 2- 
component bulge/disk model , final multi-component model) the 
2D profile plot 
Residual plot 
ID profile plot 
Model-components plot 


Again, clicking on the plots displays the full size plots (Figs. [40| ). 



Fig. 35.— Screenshot of Pipeline 4 index page. 
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S4G PIPELINE4 (101112) decomp0063.html 


http://www.oulu.fi/astronomy/ 


GAL# 633 


(prey) 

NGC1415 

(next) 


(index) 

(main) 




1-COMP SERSIC NGC1415_b 



BULGE-DISK NGC1415_bd 




FINAL NGC1415_bclbar 




Fig. 36.— Screenshot of various decomposition models for this particular galaxy. Clicking on the 
image icons opens the enlarged image. 
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INITFILE= NGC1415_bdbar.outgal 
NGC1415.phot. 1 _nonan.fits 
NGC1415_sigma2014.fits 
PSF-1.composite.fits 
CONSTRNT= none 
NGC1415.1 .finmask_nonan.fits 


FITSECT = ’[178:778,228:828]’ 
CONVBOX =’50, 50 ’ 

MAGZPT = 21.097 

INFILE: 2014-Jun-19 
PLOT: 19-Jun-2014 12:54:36.00 
s4g@fys091005 



COMP_1 = sersic 
1_XC = [478.2250] 
1_YC = [527.6920] 
1_MAG = 11.8529 
1_RE = 6.9304 
1_N = 0.9967 

1_AR = 0.4776 
1_PA = -26.7537 

COMP_2 = expdisk 
2_MAG =11.8510 
2_RS = 77.4527 
2_AR = [0.4025] 
2_PA = [-27.6820] 

COMP_3 =ferrer2 
3_MU_0 =20.2426 
3_ROUT = 91.0054 
3_ALPHA = [2.0000] 
3_BETA = [O.OOOe+00] 
3_AR = 0.4469 
3_PA = -46.0224 

Chi2/nu= 5.73102 


Fig. 37.— Final pipeline decomposition for NGC 1415. The decomposition includes bulge, disk, and 
bar components, indicated by different colors. The numbers after the labels indicate the relative 
fraction of light in each component. This 2D-profile indicates the brightness of each pixel versus 
its distance from the galaxy center. The frame also collects the names of the input data files and 
the final GALFIT decomposition parameter values. 
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Fig. 38.— Final pipeline decomposition for NGC 1415, corresponding to Fig. The upper frames 
show the observed, clipped and sky-subtracted image (left) and the corresponding masked image 
(right). The lower frames display the model image (left), and the observed-model residual image 
(gray scale covers ±1 mag). 
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NGC1415_bdbar.outgal 



sma (pix) 


Fig. 39.— One dimensional profiles corresponding to Fig. The symbols indicate the azimuthally 
averaged surface brightness as a function of semi-major axis, obtained with IRAF ellipse-routine. 
The curves indicate the azimuthally averaged profiles of the model components, using the same 
isophot es. 























Fig. 40.— Illustration of the model components corresponding to Fig. The upper left shows the 
observed image, together with superposed ellipses illustrating the various components of the final 
model: the size of the ellipse corresponds to 2 effective radii, the orientation corresponds to the 
components’ axial ratio and PA. In the upper right, the same but showing the model image. The 
lower frames are similar, except that the observed galaxy and the model have been deprojected, 
using the assumed disk orientation parameters. 
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